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ABSTRACT 


The numerical solution of differential equations of 
the form y* = f(x,y) using predictor -corrector multi- 
step methods is examined with particular emphasis on 
the stability concepts. Computational methods for de- 
termining the region of stability for single multistep 
methods and predictor-corrector pairs are expounded, 
and two subroutines have been written to compute the 
boundary of the region of stability for the single multi- 
step methods and predictor -corrector pairs. 
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stability theory of multistep methods 

By S. M. Keathley and T. J. Aird* 
Manned Spacecraft Center 

SUMMARY 


The numerical solution of differential equations of the form y' = f(x, y) using 
predictor- corrector multistep methods is examined with particular emphasis on the 
stability concepts. Computational methods for determining the region of stability for 
single multistep methods and predictor- corrector pairs are expounded. Two sub- 
routines have been written to compute the boundary of the region of stability for the 
single multistep methods and predictor-corrector pairs. 


INTRODUCTION 


The results of recent studies concerning predictor-corrector multistep methods 
for numerically solving differential equations of the form y' = f (x, y) are investigated 
with particular emphasis on the stability concepts, including justification of the theory 
for systems of differential equations. Determination of the region of stability for sin- 
gle multistep methods and predictor-corrector pairs by computational methods is pre- 
sented as outlined by Crane and Klopfenstein (ref. 1) and Krogh (ref. 2). 

A FORTRAN IV subroutine STBL1 has been written to compute the boundary of the 
region of stability for single multistep methods of either the predictor or corrector type, 
and is described along with examples of the Adams- Bashforth predictors and the 
Adams-Moulton correctors. A FORTRAN IV subroutine STBL2 has been written to 
compute the boundary of the region of stability for predictor-corrector multistep meth- 
ods for which the corrector is applied only once and is described along with examples 
of the Bashforth- Moulton pairs. 
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SYMBOLS 


f n+£ " f ( X n+f y n+c) 
h step size 

h = Ah 


W - < x n * th ) 

h * y W 

( y f a 0 

€ a = h ‘ y £ 

' t ’ y 't • 


approximate solution to y* = f (x, y) at x = x n+ ^ 
solution to y’ = f(x, y) at x = x £ 

coefficients for a multistep method 



a 

terms of order h 

DEVELOPMENT OF THEORY 


Multistep Methods 

A multistep method for solving the initial value problem y’ = f(x, y), y(x Q ) = y Q 
can be represented by the following differential -difference equation: 


G k^n+k + a k-l y n+k-l + ‘ ‘ + Vn = h (¥n + k + 'Wn+k-l + * * * + Vn) (1) 
where h is the step size, x n = x Q + nh, Y n+j = y( x n+ j)> and f n+ j = f ( x n + j’ y n + j)' 
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In order that multistep methods be useful in the numerical solution of differential equa- 
tions, it is necessary that the value of the expression 


L [y (x) ] - S, Vnrf - h £ V'n.l 


( 2 ) 




fi=0 


be small when h is small (ref. 3). By expanding the terms of Lmx)| in a Taylor 
series about x in powers of h, 


L[y(x)] = c Q y(x) + c 1 hy , (x) + . . . + c n h"^ n \x) + 


(3) 


where 


c 0 = a 0 + a l + *" +a k 


C 1 = a 1 + 2« 2 + . . . + ka fc - ^/3 Q + /3 1 + . . . + 


* lk (“l + 2 % + • • • + kd ”k) - + 


+ . . . + k 


d- 




For d = 2, 3, ... , it can be seen that 


L[y(x)] = o(h d+1 ) 


(4) 


if and only if c^ = c^ = . . . = c^ = 0; but c^ + j ^ 0. If this is the case, then d will be 

called the degree of the multistep method, whereas the integer k will be called the 
order of the multistep method. 

Example - The Adams- Bashforth k = 1, d = 1 method is 


y . - y = hf 
J n+1 n n 


(5) 
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and 



so that Cq = Cj = 0, Cg = 2* 

The corresponding Adams -Moulton k = 1, d = 1 method is 



( 7 ) 
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Stability for a Single Differential Equation 


Consider the first order equation y* = f(x, y) 


computed by using the multistep method and 


A A. 

y r y j 


and let y., y'. denote the solution 

J J 

denote the exact solution for 


j = n, . . . , n + k. Then, by substituting into the multistep formula first the computed 
values and then the exact values and subtracting, the following difference equation for 
the error is obtained. 


a k e n+k + a k-l e n+k-l 

+ . . . + a 

0 c n 

h ( /3 k c 'n+k + 1 

(The local truncation error 

L P( x n)] 

is 

assumed to be 

neglected. ) 



where 

€ r 

Y 



V 

y, j 



j = 

n, 

. . . , n + k 


Assuming that = \ is constant in the interval ^ , it follows by the mean value 
theorem that eV = Xe^, j = n, ... , n + k. Then the difference equation can be written 
as 


(“k" W k)'n+k + "- + (“0 - h,3 o)<n - 0 


( 10 ) 


where h = \h. This is a linear difference equation with constant coefficients whose 
solution may be found in the usual manner. First form the characteristic equation 
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£ k be the k roots of this equation. Then the solution to the dif- 


and let ^*2* " * * * 

ference equation is 


'n = c l«i n + c 2«2 n + -" + Vk" 


( 12 ) 


if all roots are distinct. In the case of multiple roots, let 


*j+l £j+2 " ^ 


j+m 


then 


C j+l^j+l + • • • + c j+ m ^jj+m 


is replaced by 



+ nc. 


j+2 


m- 


+ n 




where c^, ... , c^ are constants determined by the initial conditions. In either case 

(distinct or multiple roots) if | > 1 for some l, then it is clear that e will be 

large for large n; that is, the error introduced at some stage of the solution process 
will grow at each successive stage of the process until finally the solution becomes 
useless. On the other hand if | ^ | < 1 for 1 = 1, 2, ... , k, then the error e n will 

be small for large n. The above analysis leads to the following definition: Let S* be 

an arcwise connected region in the complex plane which has zero as one of its boundary 
points. Then S* will be called the region of stability for 



a k y n+k + a k-l y n+k-l + ’ ’ " + Vn 
' h (^k y ’n+k + ^k-l y ’n+k-l + ' * * + Vn) 


(13) 
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provided all roots of equation (11) have an absolute value of less than 1 for all heS*. 
If S is any other region satisfying these conditions, then S is a subset of S*. 

Ideally, the region of stability should be the left half plane. However, in prac- 
tice this is seldom the case. The Adams-Moulton k = 1, d = 2 method 


y n + l - y n = 3( y 'n + l +y 'n) <14) 

has this ideal region of stability since the characteristic equation 

('-l) 2 - M)=° 

has root z = ^ + ~ which has an absolute value of less than 1 when fte(h) < 0. This 
2 - h 

concept of stability can best be understood by considering the first order linear differ- 
ential equation y T = xy with initial condition y(0) = 1 where x is a constant. If this 
equation is substituted into the multistep formula 

Vn + k + • ' • + Vn * K (' 3 k y n + k + • • • + Vn) (15) 

and if % 2> ■ • are the roots of the characteristic equation, it would be desir- 

able to have ~ exp (hA) and (jg = • • • = 4^ = 0. Then < 1 when Re(A) < 0, 
and |£i| > 1 when Re(A) > 0. That is to say, the multistep method has the left half 
plane as its region of stability. 

Stability for a System of Differential Equations 

Consider the system of differential equations Y’ = F(x, Y) where Y = (y^, y 2 , • • •> 
y n ^ and F = (f j, i 2> • . . , f n ). Let Y^x^, Y’^c. ^ denote the solution computed by using 
the multistep method and Y^x^, denote the exact solution for j = n, . . . , n + k. 
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Then by substituting into the multistep formula first the computed values and then the 
exact values and subtracting, the following difference equation for the error is obtained 


V( X n + k) * “k-l« fn + k-l) + ' ' ' + “o'( X n) * h [V( x n + k) ♦ « k -l'’( X n4k-l) 

+ . . . + S 0 ''( X n)] < 16 > 

(The local truncation error L Y ^x A is assumed to be constant and is therefore 
neglected. ) 

where f j) = Y ( x j) ' Y {*j) 

e '( X i) = Y '( X i) ' Y ( X i) 


j = n, ... , n + k 
3f. 

Let A = fa. ^ where a. the terms f. and y. being the ith and jth components 

V ij) 1 3 

of the vectors F and Y, respectively. Assuming that a„ is constant, it follows by 

the mean value theorem that e' = Af. Then the difference equation can be written as 


“k'( x n+k) + “k-l'( x n+k-l) + • • * 


+ “o f ( x n + k-l) 

- + \-l'( x n + k-l) + • • • + V( x ,| 


(17) 


Since every square matrix is similar to its Jordan canonical form, there exists a non- 
singular matrix P such that PAP ^ = J(A) where 
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is the Jordan canonical form of A and 



where X., 1 = 1, ... , m are the eigenvalues of A and need not all be distinct. Then 

1 _i 

let E = Pe, solve for e = P E, substitute into the difference equation, and multiply 
both sides by P to obtain 

“k E ( x „ + k) + • • • + °0 E (*n) ■ hJ < A >fk E ( x „ + k) + • ' • + V 2 ( x n)] < 20 > 

The solution to each equation of this system of difference equations can be written in 
terms of powers of the roots of the characteristic equations 


( a k " X i h/3 k) zk + • • • + (« 0 " X i h ^ 0 ) = 0; 1 = ■ • • 2 


For a more rigorous treatment see Lea (ref. 4). 

_ Therefore, with a system of differential equations the stability is related to a 
set hj = XJi where X. are the eigenvalues of the Jacobian matrix of the derivative 

functions. These eigenvalues may be complex numbers so that the stability must be 
considered in terms of complex values of h. 

Let i-g, H = 1, ... , k, be the roots of the characteristic equation (21). Then if 

Jl^j < 1, S. = 1, ... , k, the vector E^x n J will be small for large n and therefore 
the error vector e(x n ) = P -1 E(x n ) will also be small. 

The stability definition for a single equation will suffice for systems of equations. 
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Finding the Region of Stability 

The problem may be restated as follows: For what values of h in the given 

characteristic equation 



( 22 ) 


do the roots have an absolute value of less than 1. The direct method of finding the 
region of stability would be to vary h at each step, examining the roots to determine 
whether h was in the region. As k increases, the calculations involved soon become 
prohibitive. If the problem is looked at from a slightly different viewpoint, the answer 
is forthcoming. Because of the existing continuity between the roots and the coeffi- 
cients of the characteristic equation, it is only necessary to calculate the values of h 
for which the equation has roots of an absolute value equal to 1. For this purpose it 
is necessary to rewrite the characteristic equation as 


k 

1=0 1=0 

which is linear in h. Thus 



h = 


i_A z 

a=o 

k t 
£=0 


(24) 


To find those values of h for which the characteristic equation has roots of an abso- 

i0 

lute value equal to 1, let z = e . Then 


h(6) 


H=0 


He 


Ev 

1=0 


iSL8 


(25) 
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By plotting h(e) for 0 < 0 < n and the conjugate, the region of stability is easily 
found. In any arcwise connected region, if the method is stable at one point, then it fol- 
lows that it must be stable in the entire region; that is, h must cross the boundary 
h(e) in order to go from a stable region to an unstable region. 

Example - The Adams -Bashforth k = 1, d = 1 method 


y n + l 


(26) 


shown in sketch a has the characteristic equation 



Sketch a 


The Adams-Moulton k = 1, d 


1 method 


z - (1 + h) = 0 


Therefore, 


h<6) = e ie - 1 


( 27 ) 


(28) 


y n+ l" y n = hy, n + l 


(29) 


shown in sketch b has the characteristic equation 



Sketch b 


11 



A FORTRAN IV subroutine STBL1 has been written to compute the boundary of the re- 
gion of stability for single multistep methods. The subroutine is described in appen- 
dix A along with examples of the Adams- Bashforth predictors and the Adams-Moulton 
correctors. 


Stability of Predictor- Corrector Pairs 
The following two multistep methods define a predictor-corrector pair 

L p[ y ( X -)] = a k y n+k + a k-lVk-l + ' • • + Vn 

~ h (^k-l^'n+k-l ' f ' ' ‘ +Vn) < 32 > 

L c[ y ( x n)] - v kVk + >k-l y n + k-l + • • • + ’Vn 

- h (V.* + 5 k-l y 'n+k-l + • • • V’n) (33 > 


The predictor is used first to give a value for y n+ j J > which is then used to 

evaluate y’^^, which in turn is used by the corrector to compute a new value for 

y n+ ^. This process may be iterated; however, usually the degrees of and L c 

are equal, and it is assumed that convergence is reached after one application of the 
corrector. The difference between y n+ ^ of the predictor and y n+ ^ of the corrector 

is used as an estimate of the local truncation error (ref. 5). In particular, 


J p[ y ( X n)] = 


= c 


p, d+1 


h dtl y (d + l)^ + 0 ( h <d + 2)) 


(34) 


and 




C c, d+l h 


d+1 


y^) ( x n> + 0(h( d+2 )) 


(35) 


It follows that 


z, n+k " y p, n+k- ( C c, d+1 ' C p, d+l) h d+ V d+1) (x n ) 


(36) 
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or 


'c, d+1 


-h 


(37) 


The advantage of such an error estimate makes it worthwhile to have the degree of the 
predictor equal to the degree of the corrector. 

The characteristic equation of the pair can be derived (ref. 6) by letting y' = Ay, 
h = Ah. Then from the predictor 



- hS k-l) y n^-l + " 



(38) 


By substituting this value into the corrector (assuming that a k = 1) 


Vn+k + y k-l y n+k-l + ’ * ’ + y 0 y n 

= h {-«k[(“k-l -^k-l) y „ + k-l + • • • + (“o -^oK] 

+ «k-l y n + k-l + • • • + Vnj (39) 

and then the characteristic equation is 

( r k) zk + [ y k-i +lifi k(“k-i -H-i) -"^k-ij 11 ' 1 

+ • • • + [j0 - h8 0 ) -Ti« 0 J = 0 (40) 
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Again it is desirable to rewrite this equation in powers of h and compute the boundary 
of the region of stability by computing the values of h for which the characteristic 
equation has roots of an absolute value equal to 1. The equation in powers of h is 


Let 



(41) 


k-1 




e ife 


c = 


-L 

1=0 


n 


me 

e 


Then 

Hi ( 0) = -B^vgTiAC ( 42 ) 

and 


h 2 te> = - B - V f A - 4AC (43) 

By plotting hj(e) and hgte) for 0< 0 < n and the conjugates, the region of stability 

for the predictor-corrector is easily found in exactly the manner described for single 
formulas. 
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Example - The Adams- Bashforth k = 1, d = 1 predictor and the Adams-Moulton 
k = 1, d = 1 corrector are 


y n+l - y n = hy 'n < 44 > 


and 


^n+1 " ^n ^ n+1 


The joint characteristic equation is 


(h) 2 + h - e 10 + 1 = 0 (46) 


shown in sketch c and 



Sketch c 


A FORTRAN IV subroutine STBL2 has been written to compute the boundary of 
the region of stability for predictor-corrector multistep methods. The subroutine is 
described in appendix B along with examples of the Bashforth- Moulton pairs. 

The generalization of this method to finding the region of stability of predictor- 
corrector pairs where the corrector is iterated m times _is apparent. The charac- 
teristic equation will be a polynomial of degree m + 1 in h. The details will not be 
discussed here. 
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CONCLUDING REMARKS 


Algorithms have been developed for calculating the region of stability for single 
multistep methods and for predictor-corrector pairs. The results are easily general- 
ized to include other versions of the predictor-corrector methods. 

It should be noted that, in both cases considered, an explicit formula is given for 
the boundary of the region of stability. This formula may be useful in developing a 
technique to generate multistep methods with an increased region of stability. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, January 19, 1967 
039-00-00-00-72 
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APPENDIX A 


FORTRAN IV SUBROUTINE STBLl FOR SINGLE MULTISTEP METHODS 


The subroutine STBLl computes the boundary of the region of stability for a sin- 
gle multistep method of either the predictor or corrector type. 

In order to use the subroutine STBLl, the following calling sequence is necessary: 

CALL STBLl (AA, BB, K, N, TITLE) 

Where the multistep method is 

Vn+k * Vl y n+k-l + • ' ' * “Qyn = h (Vn + k + 'Wn+k-1 + • • • + ¥i) 


and 


AA is a double precision array with 25 locations such that AA(1) = 

AA(2) = , AA(k + 1) = a Q . 

BB is a double precision array with 25 locations such that BB(1) = 

BB(2) = ) 3 k _ 1 , ... , BB(k + 1) = / 3 q . 

K is an integer giving the order of the method. 

N is an integer giving the number of divisions of the interval (0, n) to be used 
in computing the boundary of the region of stability; 90 divisions are usu- 
ally sufficient. 

TITLE is an array of 5 locations containing BCD information to be used as a title 
for the printed and plotted output. TITLE may be passed to the subroutine 
as 30H (any 30 BCD characters). 

STBLl will compute the boundary of the region of stability and print the results 
on the line printer if 10 = 6 or on the 4020 output if 10 = 17; it will also plot the results 
on the 4020 plotter. 

Table AI is a listing of the FORTRAN IV subroutine STBLl. Subroutine STBLl 
has been used to calculate the regions of stability for the Adams-Bashforth predictors 
and the Adams-Moulton correctors of degree 1 to 4. Table AH gives the coefficients 
of the Bashforth methods and the numerical results of STBLl, while figure Al displays 
the results graphically. Table Am gives the coefficients of the Moulton methods and the 
numerical results of STBLl, while figure A2 displays the results graphically. The 
coefficients for both of these methods are available from NASA Manned Spacecraft 
Center, Houston, Texas. 
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TABLE AI. - LISTING OF FORTRAN IV SUBROUTINE STBL1 FOR 
SINGLE MULTISTEP METHODS 


DIMENSION AAl2b)f bBi?5) 

DOUBLE PRECISION Ah* RB 
aAU) _= l.uDO 
AA ( 2 ) = -1 .OUO 
UO. Iu2b o = 3»25 
lU2b AAIJ) = U.UDu 
tiB < 1 ) = U.UDU 
HR ( 2 ) = i.uDu 

KSTEP = 1 _ _ 

CALL STRLlt AA rdBfKSTe r» r luOr 
1 3UHADAMS b A SHE ORTH t K=l» U=1 

bRU) = U.UDu 

Rb(2) = U. ISOOOOUOUOOOOOUODI 

Rb<3> = -u.souoooooonoooououn 
kSTEP = 2 

CALL STBLllAA.bR, KbTtP*lljOf 
1 3UHADAMS bASHEOKTH t K=2# 0=2 

BR ( 1 ) = U.UDU 

Rb<2> = u. iqi6b6bfeb6b6b6b7L)l 
RB(3) = -u .133333333333333301 
Rb(4) = U . 4J b6fo6b6bf>b6b6b7U0 
rSTEH '= 3 

CALL STBL1 i AA.BP.KSTtP* 100. 

1 3UHADAMS bASHEOKTH, K=3.* 0=3 

dEU 1 ) = U.UDu 

Rb (2) = U • 229 Ib6bbo666b6670 1 

Rb (3 ) = -U.24bB333333333333Ul 

RB(4) = U.'l'S41b6bbb6b666b7L)1 

BB(S) = -U.37bO(jOuOOOUOoOuOL)n 
kSIEP = 4 

call STRLlt A a.br, KsTlP* luO, 

1 3UHADAMS BASHEOKTH. k=4, 0=4 

ciBtlJ = i.UDu 

bR ( 2 ) = U.UDU 

rSIEP = 1 

CALL STBLUAA, BB, KSTlP* luO, 

1 3UHADAMS MOULTON K=1 » D=i 

Rb ( 1 ) = 0*500 
Pb ( 2 > = 0.500 
kSTEH = 1 

CALL STBL1 tAA.BR, KSTtP* luO , 

1 3UHADAMS MOULTON K = 1 , 0=:? 

Rti(l) = 0 • 4 1 6b6b6b6b566bB7D— 0 0 
Hb(2) = 0 * 6b6b6b6b6b666b67D— n U 
_ Rb ( d > =- 0 .8333333333333333D-Q1 
KSfEH = d 

CALL STBLl(AA,BB,KsTtP.lUO. 

1 3UHADAMS MOULTON K = 2, D =3 

Pb ( 1 ) = U.373D0 

Rb(2) = U • 7Q16b6b6bBb66bb7U-U0 

Bb<3> = “U • 2 0b33333 33333 3 33U- 00 

Rb(4) — U . h 1 b6b6b6b65bb6b7u— U 1 

rvSfEH = 3 

CALL STBLl (AA,BR,KsTtP.lt)0. 

1 30HADA^b MOULTON * = 3, 0=4 

tNU 


) 


) 


) 


) 


) 


) 
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TABLE AI. - LISTING OF FORTRAN IV SUBROUTINE STBLl FOR 


SINGLE MULTISTEP METHODS - Concluded 


SUBROUTINE STRLKAm* Rbf Kf N* TITLE) 

U I MENS ION mA (? b) * bH(Pb), hI(rOO)* H^fbOo)* TTTLF<5> 

UOUBLE PkEcISION Aa* RB* XI* X2* Yl. Y** OblN* DcOS* J)EL* THETA 
UOUHLE PRECISION DL * PI * 02* DL) 

DIMENSION bCux(12>* t»CDY (12) 

0 1 MENS ION im5U0) 

UATA BCDX/12*6H / 

L) A 1 A RCn>/i2*6H / 

10=6 
10=1 1 
M=K + 1 
1C=0 
NC=0 

L>EL=3. 1415y2t,5^5ti979bP0/N 

THLTA=0-UDU _ 

1U2U IF (NC .t»T. hi) GO TO 1030 
IHC.TA = NC»DtL 

Xl=O.ODO 
X2=0.0U0 
Y 1 =0 ■ 0L>0 

Y2=0 - nuo 

UO 1U1U 1 = 1 * K 1 

L=N-i+l 

UL=L 

Xl=Xl+AAi I J*UC0S(DL*1 HETa) 

X?=X2+AA l I > ♦US1 N(DL*IHETa) 

Yl=Tl + bBU)*UCUS(DL^THLTA) „ . 

Y?=Y2+bBl I ) *USIN (DL*1 HLTa ) 

1010 continue 

UD=Yl**2+Y2**2 

IF 10U .LI. 1.0U-b> Gu TO 1040 

ul=< xl*Yl+A2*Y2) /DU 

U2=<X2*Y1-a1*Y2>/DU 

IC=IC+i 

HlUO=Dl 

H2( IC)=D^ 

IHIIC) = ThETA 
1040 NC=NC+1 

GO LQ_ 1CU0 

1030 CONTINUE 

6001 h0RPATUhl»40Xf4K — " rSA b r4H *///) 

rtPITt < 10* 6u0l ) (TIlLt- ( J) f J=1 rS) 

UO 1060 o=l#Kl 
L=K-U+1 

1060 WRITE ( 10p 6UO^) L* AA( J) .#L, Bb< J) 

6U02 FOKMATtlN * lbX *6H ALPHA ( * I2*4H) = *D21.16.?0X* 

1 6h bt |‘A ( * 12 * 4H) = f L>2 1 • 16 ) 

6U0b P OKM AT {// / i 

rtRlTE< 10* 6u0b) 

6003 POKMATUh » l3X*E10.4,sX* 1HI *L1 0.4*3H p F 1 0 . 4 p I H ) r 1 OX , 

_ 1 .. .. f!0,._4»5X*JLH( *tlU.A*3H * El0.4»lH)) 

6U0t) POKPATUh * lbX*5HThEl A» lMX»4HMBAKr23X»bHTHETA» IPX f 4HHRAP/) 

WRITE ( 10 » 6oO b ) 

1C2=( IC+1 >/2 
UO 1U7U 1=1*1C2 
J= 1+ IC/2 

YIP 6 OO 3 ) < 1H_1_I ) pHI ( I) *H£( I > * TH( S U *H1 ( J)_*H2< J) ) 

1070 CONTINUE 

UO 1U8U i=i*6 
10P0 nCuxn ) = I 1 1 LE ( 1 ) 

CALL QUlt\ML(-lp-3.0'3.0*-3.0f3.nplH r RCDX * UCDY , IC * Hi * H2 ) 

UO 10R0 i=l*lC 
luou .h.pl I ) =-N2U ) 

CALL OUlNML(up-3.n#3.0* -3. 0*3.0* 1H p uCUX p RCDY * IC p HI * H2 ) 

HE I UKN 
ENU 
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TABLE AH. - LINE PRINTOUTS OF ADAMS -BASHFORTH PREDICTORS 

(a) k = 1, d = 1 


ALPHA ( 1) = 
ALPHA { 0 ) - 


THETA 

• 0000 
•3142-01 
•6283-01 
.9425-01 
.1257-00 
•1571-00 
,1885-00 
•2199-00 
•2513-00 
.2827-00 
•3142-00 
•3456-00 
•3770-00 
•4084-00 
,4398-00 
.4^712-00 
•5027-00 
,5341-00 
•5655-00 
•5969-00 
•6283-00 
•6597-00 
.6912-00 
•7226-00 
•7540-00 
•7854-00 
•8168-00 
.8482-00 
,8796-00 
•9111-00 
.9425-00 
•9739-00 
,1005+01 
•1037+01 
.1068+01 
• 1100+01 
•1131+01 
•1162+01 
.1194+01 
•1225+01 
•1257+01 
•1288+01 
•1319+01 
,1351+01 
•1382+01 
,1414+01 
.1445+01 
,1477+01 
.1508+01 
•1539+01 
•1571+01 


ADAMS BASHFORTH f K=l, D=1 


• 1000000000000000+01 

-. 1000000000000000+01 


BETA ( 1) = .0000000000000000 

BETA ( 0) = .1000000000000000+01 


HBAR THETA HBAR 



• 0000 ' 

-.0000 ) 

,1571+01 

( -.1000+01 f 

.1000+01) 


-•4934-03 f 

•3141-01) 

.1602+01 

( -.1031+01 r 

.9995-00) 


-.1973-02 f 

•6279-01) 

.1634+01 

( -.1063+01 f 

.9980-00) 


-.4438-02 f 

.9411-01) 

•1665+01 

( -.1094+01 f 

,9956-00) 


-•7885-02 f 

•1253-00) 

,1696+01 

( -.1125+01 * 

.9921-00) 


-.1231-01 f 

.1564-00) 

.1728+01 

( -.1156+01 » 

.9877-00) 


-.1771-01 f 

.1874-00) 

.1759+01 

( -.1187+01 * 

.9823-00) 


-.2408-01 f 

.2181-00) 

,1791+01 

( -.1218+01 # 

.9759-00) 


-.3142-01 f 

.2487-001 

.1822+01 

( -.1249+01 f 

.9686-00) 


-.3971-01 r 

.2790-00) 

.1854+01 

( -.1279+01 * 

.9603-00) 


-.4894-01 f 

.3090-00) 

•1885+01 

( -.1309+01 f 

.9511-00) 


-.5912-01 f 

.3387-00) 

.1916+01 

< -.1339+01 r 

.9409-00) 


-.7022-01 f 

.3681-00) 

.1948+01 

< -.1368+01 r 

.9298-00) 


-.8225-01 * 

.3971-00) 

.1979+01 

( -.1397+01 » 

.9178-00) 


-•9517-01 » 

.4258-00) 

,2011+01 

( -.1426+01 > 

.9048-00) 


-.1090+00 r 

.4540-00) 

,2042+01 

( -.1454+01 f 

.8910-00) 


-•1237+00 f 

.4818-00) 

•2073+01 

( -.1482+01 > 

.8763-00) 


-.1393-00 f 

.5090-00) 

.2105+01 

( -.1509+01 > 

.8607-00) 


-.1557-00 f 

.5358-00) 

.2136+01 

( -.1536+01 » 

.8443-00) 


-.1729-00 » 

.5621-00) 

.2168+01 

( -.1562+01 f 

.8271-00) 


-.1910-00 f 

.5878-00) 

•2199+01 

( -.1588+01 » 

.8090-00) 


-•2098-00 f 

.6129-00) 

.2231+01 

( -.1613+01 » 

.7902-00) 


-.2295-00 r 

.6374-00) 

.2262+01 

( -.1637+01 » 

.7705-00) 


-.2499-00 f 

.6613-00) 

•2293+01 

( -.1661+01 f 

.7501-00) 


-•2710-00 r 

.6845-00) 

.2325+01 

< -.1685+01 f 

.7290-00) 


-.2929-00 f 

.7071-00) 

.2356+01 

< -.1707+01 r 

.7071-00) 


-.3155-00 * 

•7290-00) 

.2388+01 

{ -.1729+01 f 

.6845-00) 


-.3387-00 r 

•7501-00) 

•2419+01 

( -.1750+01 f 

.6613-00) 


-•3626-00 r 

.7705-00) 

.2450+01 

( -.1771+01 f 

.6374-00) 


-.3871-00 f 

.7902-00) 

.2482+01 

( -.1790+01 f 

.6129-00) 


-•4122-00 r 

.8090-00) 

.2513+01 

( -.1809+01 f 

.5878-00) 


-.4379-00 f 

.8271-00) 

•2545+01 

( -.1827+01 f 

.5621-00) 


-.4642-00 f 

.8443-00) 

.2576+01 

< -.1844+01 f 

.5358-00) 


-•4910-00 t 

.8607-00) 

.2608+01 

( -.1861+01 f 

.5090-00) 


-.5182-00 f 

,8763-00) 

,2639+01 

( -.1876+01 f 

.4818-00) 


-.5460-00 f 

.8910-00) 

,2670+01 

( -.1891+01 f 

.4540-00) 


-.5742-00 r 

•9048-00) 

.2702+01 

< -.1905+01 f 

.4258-00) 


-.6029-00 f 

.9178-00) 

•2733+01 

( -.1918+01 f 

.3971-00) 


-.6319-00 f 

.9298-00) 

•2765+01 

( -.1930+01 f 

.3681-00) 


-.6613-00 r 

.9409-00) 

.2796+01 

( -.1941+01 f 

.3387-00) 


-.6910-00 » 

.9511-00) 

.2827+01 

( -.1951+01 f 

,3090-00) 


-.7210-00 * 

.9603-00) 

.2859+01 

( -.1960+01 f 

.2790-00) 


-.7513-00 f 

.9686-00) 

.2890+01 

( -.1969+01 f 

.2487-00) 


-.7819-00 f 

.9759-00) 

•2922+01 

( -.1976+01 f 

.2181-00) 


-•8126-00 f 

•9823-00) 

•2953+01 

( -.1982+01 f 

.1874-00) 


-.8436-00 t 

•9877-00) 

.2985+01 

( -.1988+01 f 

,1564-00) 


-.8747-00 f 

.9921-00) 

.3016+01 

( -.1992+01 f 

,1253-00) 


-.9059-00 f 

•9956-00) 

.3047+01 

< -.1996+01 f 

.9411-01) 


-•9372-00 f 

.9980-00) 

•3079+01 

( -.1998+01 f 

.6279-01) 


-.9686-00 f 

•9995-00) 

.3110+01 

( -.2000+01 f 

,3141-01) 


-.1000+01 f 

.1000+01) 

.3142+01 

( -.2000+01 f 

.9992-15) 



TABLE AH. - LINE PRINTOUTS OF ADAMS -BASH FORTH PREDICTORS - Continued 

(b) k = 2, d = 2 


ADAMS BASHFORTHf K=2f 0=2 


ALPHA ( 2) = .1000000000000000+01 
ALPHA ( 1) = -.1000000000000000+01 
ALPHA ( 0) = .0000000000000000 


BETA< 2) = .0000000000000000 

BETA< 1) = .1500000000000000+01 

BETA< 0) = -.5000000000000000+00 


THETA 

HBAR 

THETA 

HBAR 


,0000 

( ,0000 f 

-.0000 ) 

.1571+01 

( -.4000-00 

F 

.8000-00) 

.3142-01 

( -.2433-06 r 

.3140-01) 

.1602+01 

( -.4177-00 

f 

,7971-00) 

,6283-01 

< -.3882-05 f 

.6273-01) 

.1634+01 

( -.4354-00 

F 

.7936-00) 

.9425-01 

( -.1957-04 f 

.9390-01) 

.1665+01 

( -.4532-00 

F 

.7894-00) 

,1257-00 

( -.6145-04 f 

.1248+00) 

.1696+01 

( -.4711-00 

F 

.7844-00) 

.1571-00 

< -.1488-03 f 

.1555-00) 

.1728+01 

( -.4890-00 

F 

.7789-00) 

,1885-00 

( -.3056-03 f 

.1858-00) 

.1759+01 

( -.5070-00 

F 

.7726-00) 

,2199-00 

( -.5598-03 f 

,2156-00) 

.1791+01 

( -.5249-00 

F 

.7657-00) 

,2513-00 

( -.9426-03 f 

,2450-00) 

.1822+01 

< -.5427-00 

F 

.7581-00) 

.2827-00 

( -.1488-02 f 

.2738-00) 

.1854+01 

< -.5605-00 

F 

.7499-00) 

,3142-00 

t -.2232-02 f 

.3020-00) 

.1885+01 

< -.5782-00 

F 

.7410-00) 

,3456-00 

( -.3210-02 f 

,3295-00) 

,1916+01 

( -.5958-00 

F 

.7315-00) 

,3770-00 

( -.4461-02 f 

,3564-00) 

.1948+01 

( -.6133-00 

F 

.7214-00) 

,4084-00 

( -.6021-02 f 

.3826-00) 

,1979+01 

( -.6306-00 

F 

.7107-00) 

.4398-00 

( -.7926-02 f 

.4080-00) 

.2011+01 

C -.6477-00 

F 

.6993-00) 

.4712-00 

( -.1021-01 f 

.4327-00) 

.2042+01 

( -.6646-00 

F 

.6874-00) 

,5027-00 

( -.1291-01 f 

,4566-00) 

,2073+01 

( -.6813-00 

F 

.6748-00) 

.5341-00 

< -.1604-01 f 

.4797-00) 

.2105+01 

( -.6978-00 

F 

,6617-00) 

.5655-00 

( -.1965-01 f 

.5020-00) 

,2136+01 

( -.7140-00 

F 

.6481-00) 

.5969-00 

C -.2374-01 f 

,5235-00) 

.2168+01 

{ -.7299-00 

F 

,6339-00) 

,6283-00 

( -.2835-01 f 

.5442-00) 

.2199+01 

( -.7455-00 

F 

.6191-00) 

.6597-00 

( -.3349-01 f 

.5640-00) 

.2231+01 

< -.7608-00 

F 

.6038-00) 

,6912-00 

( -.3918-01 f 

.,5830-00) 

.2262+01 

( -.7758-00 

F 

.5880-00) 

.7226-00 

( -.4542-01 f 

,6012-00) 

,2293+01 

( -.7904-00 

F 

.5717-00) 

,7540-00 

( -.5223-01 f 

.6186-00) 

.2325+01 

( -.8046-00 

F 

,5549-00) 

,7854-00 

( -.5960-01 f 

.6352-00) 

,2356+01 

( -.8184-00 

F 

.5376-00) 

.8168-00 

C -.6755-01 f 

.6509-00) 

,2388+01 

< -.8319-00 

F 

.5199-00) 

,8482-00 

( -.7607-01 f 

.6659-00) 

,2419+01 

( -.8449-00 

F 

.5017-00) 

,8796-00 

t -.8515-01 f 

,6800-00) 

,2450+01 

( -.8575-00 

F 

.4831-00) 

.9111-00 

( -.9480-01 f 

.6934-00) 

.2482+01 

< -.8696-00 

F 

.4640-00) 

,9425-00 

( -.1050+00 f 

.7060-00) 

.2513+01 

( -.8812-00 

F 

.4446-00) 

.9739-00 

t -.1157+00 f 

,7178-00) 

,2545+01 

( -.8924-00 

F 

.4248-00) 

.1005+01 

( -.1270-00 f 

.7288-00) 

,2576+01 

( -.9031-00 

F 

.4046-00) 

.1037+01 

{ -.1388-00 f 

.7391-00) 

.2608+01 

< -.9133-00 

F 

.3841-00) 

,1068+01 

( -.1511-00 f 

,7485-00) 

.2639+01 

( -.9229-00 

F 

.3633-00) 

.1100+01 

( -.1639-00 f 

.7573-00) 

.2670+01 

< -.9321-00 

F 

.3421-00) 

,1131+01 

( -.1771-00 f 

.7653-00) 

,2702+01 

( -.9407-00 

F 

.3206-00) 

.1162+01 

( -.1908-00 f 

.7725-00) 

.2733+01 

( -.9487-00 

F 

.2989-00) 

.1194+01 

( -.2050-00 f 

,7790-00) 

.2765+01 

( -.9562-00 

F 

.2769-00) 

.1225+01 

( -.2195-00 f 

.7847-00) 

.2796+01 

( -.9631-00 

F 

.2547-00) 

.1257+01 

( -.2345-00 f 

.7897-00) 

,2827+01 

< -.9694-00 

F 

.2322-00) 

.1288+01 

( -.2497-00 f 

.7940-00) 

.2859+01 

( -.9752-00 

F 

.2096-00) 

.1319+01 

( -.2654-00 f 

.7975-00) 

.2890+01 

( -.9804-00 

F 

•1868-00) 

,1351+01 

( -.2813-00 f 

.8003-00) 

.2922+01 

( -.9850-00 

F 

.1638-00) 

.1382+01 

( -.2976-00 f 

.8024-00) 

,2953+01 

( -.9889-00 

F 

•1406-00) 

.1414+01 

( -.3141-00 f 

.8038-00) 

.2985+01 

< -.9923-00 

F 

.1174+00) 

,1445+01 

( -.3309-00 f 

.8044-00) 

.3016+01 

( -.9951-00 

F 

.9403-01) 

,1477+01 

( -.3479-00 f 

.8044-00) 

.3047+01 

C -.9972-00 

F 

.7059-01) 

,1508+01 

( -.3651-00 f 

.8036-00) 

,3079+01 

( -.9988-00 

F 

.4710-01) 

.1539+01 

( -.3825-00 f 

.8022-00) 

.3110+01 

( -.9997-00 

F 

.2356-01) 

,1571+01 

( -.4000-00 f 

,8000-00) 

,3142+01 

( -.1000+01 

F 

.7494-15) 



TABLE AH. - LINE PRINTOUTS OF ADAMS -BASH FORTH PREDICTORS - Continued 


(c) k = 3, d = 3 





ADAMS BASHFORTH r 

K=3 1 D=3 

— 


ALPHA ( 

3) 

_ 

.1000000000000000+01 

BETAt 

3) = 

.0000000000000000 

ALPHA { 

2) 


-.1000000000000000+01 

BETA ( 

2) = 

. 1916666666666666+01 

ALPHA { 

1) 

- 

.0000000000000000 

BETA ( 

1) = 

-.1333333333333332+01 

ALPHA ( 

0) 

= 

.0000000000000000 

BETAt 

0) = 

.4166666666666666+00 


THETA 

HBAR 

THETA 

HBAR 

.0000 

t .0000 f 

-.0000 ) 

. 1571+01 

( -. 4138-01 f 

. 7034 - 00 ) 

. 3142-01 

t . 3651-06 f 

. 3142 - 01 ) 

. 1602+01 

( -. 5513-01 f 

. 6965 - 00 ) 

, 6283-01 

( . 5830-05 f 

. 6283 - 01 ) 

, 1634+01 

{ -. 6916-01 f 

. 6891 - 00 ) 

, 9425-01 

t , 2942-04 f 

. 9424 - 01 ) 

, 1665+01 

( -. 8343-01 f 

. 6815 - 00 ) 

. 1257-00 

( . 9255-04 » 

. 1257 - 00 ) 

. 1696+01 

( - . 9790 -oi f 

. 6736 - 00 ) 

. 1571-00 

( , 2246-03 t 

. 1570 - 00 ) 

. 1728+01 

( -. 1125+00 f 

. 6653 - 00 ) 

. 1885-00 

( . 4623-03 f 

. 1884 - 00 ) 

. 1759+01 

( -. 1273-00 r 

. 6568 - 00 ) 

. 2199-00 

( , 8486-03 * 

. 2197 - 00 ) 

. 1791+01 

( -. 1421-00 f 

. 6479 - 00 ) 

. 2513-00 

l . 1432-02 f 

. 2509 - 00 ) 

. 1822+01 

( -. 1570-00 f 

. 6387 - 00 ) 

. 2827-00 

( , 2264-02 p 

. 2820 - 00 ) 

, 1854+01 

( -. 1719-00 f 

. 6292 - 00 ) 

. 3142-00 

( , 3400-02 t 

. 3129 - 00 ) 

. 1885+01 

( -. 1867-00 f 

. 6194 - 00 ) 

, 3456-00 

( . 4892-02 f 

. 3436 - 00 ) 

. 1916+01 

( -. 2016-00 f 

. 6092 - 00 ) 

, 3770-00 

( . 6792-02 f 

. 3739 - 00 ) 

. 1948+01 

( -. 2163-00 f 

. 5988 - 00 ) 

. 4084-00 

( . 9147-02 f 

. 4039 - 00 ) 

. 1979+01 

( -. 2310-00 f 

. 5879 - 00 ) 

, 4398-00 

( . 1199-01 f 

. 4333 - 00 ) 

,2011+01 

( -. 2455-00 f 

. 5768 - 00 ) 

. 4712-00 

( , 1536-01 f 

. 4621 - 00 ) 

. 2042+01 

< -. 2599-00 f 

. 5654 - 00 ) 

. 5027-00 

( . 1925-01 f 

. 4901 - 00 ) 

, 2073+01 

( -. 2741-00 f 

. 5536 - 00 ) 

. 5341-00 

( . 2367-01 f 

. 5173 - 00 ) 

, 2105+01 

( -. 2881-00 f 

. 5415 - 00 ) 

. 5655-00 

( . 2859-01 f 

. 5434 - 00 ) 

. 2136+01 

t -. 3019-00 f 

. 5290 - 00 ) 

. 5969-00 

( . 3397-01 f 

. 5684 - 00 ) 

. 2168+01 

< -. 3155-00 f 

. 5163 - 00 ) 

. 6283-00 

t . 3974-01 f 

. 5921 - 00 ) 

. 2199+01 

{ -. 3288-00 f 

. 5032 - 00 ) 

. 6597-00 

( . 4581-01 f 

. 6144 - 00 ) 

. 2231+01 

t -. 3419-00 f 

. 4898 - 00 ) 

. 6912-00 

( , 5208-01 f 

. 6352 - 00 ) 

.2262+01 

( -. 3546-00 f 

. 4761 - 00 ) 

. 7226-00 

( , 5844-01 f 

. 6544 - 00 ) 

, 2293+01 

( -. 3671-00 f 

. 4621 - 00 ) 

. 7540-00 

( , 6474-01 f 

. 6720 - 00 ) 

. 2325+01 

( -. 3792-00 f 

, 4478 - 00 ) 

, 7854-00 

( . 7085-01 f 

, 6878 - 00 ) 

. 2356+01 

( -. 3910-00 f 

. 4332 - 00 ) 

, 8168-00 

( . 7662-01 f 

. 7020 - 00 ) 

. 2388+01 

( -. 4025-00 f 

. 4183 - 00 ) 

. 8482-00 

( . 8192-01 f 

. 7144 - 00 ) 

. 2419+01 

( -. 4136-00 f 

. 4032 - 00 ) 

, 8796-00 

( . 8661-01 f 

. 7252 - 00 ) 

. 2450+01 

( -. 4243-00 f 

. 3877 - 00 ) 

. 9111-00 

( , 9058-01 f 

, 7343 - 00 ) 

. 2482+01 

t -. 4346-00 f 

. 3720 - 00 ) 

. 9425-00 

( . 9371-01 f 

, 7420 - 00 ) 

. 2513+01 

( -. 4446-00 f 

. 3561 - 00 ) 

. 9739-00 

t . 9593-01 f 

. 7481 - 00 ) 

. 2545+01 

( -. 4541-00 f 

. 3399 - 00 ) 

, 1005+01 

t . 9715-01 f 

. 7529 - 00 ) 

, 2576+01 

( -. 4632-00 f 

. 3234 - 00 ) 

, 1037+01 

< . 9734-01 f 

. 7564 - 00 ) 

. 2608+01 

t -. 4718-00 f 

. 3068 - 00 ) 

. 1068+01 

t , 9646-01 f 

, 7588 - 00 ) 

. 2639+01 

t -. 4800-00 f 

. 2899 - 00 ) 

.1100+01 

( , 9450-01 f 

. 7601 - 00 ) 

, 2670+01 

( -. 4878-00 f 

. 2728 - 00 ) 

. 1131+01 

t , 9145-01 f 

. 7605 - 00 ) 

, 2702+01 

( -. 4951-00 f 

. 2555 - 00 ) 

. 1162+01 

( . 8735-01 / 

. 7600 - 00 ) 

. 2733+01 

( -. 5019-00 f 

. 2380 - 00 ) 

, 1194+01 

( . 8222-01 f 

, 7587 - 00 ) 

. 2765+01 

( -. 5083-00 f 

. 2204 - 00 ) 

. 1225+01 

t . 7609-01 f 

. 7567 - 00 ) 

. 2796+01 

( -. 5142-00 f 

. 2026 - 00 ) 

, 1257+01 

( . 6901-01 f 

, 7541 - 00 ) 

. 2827+01 

{ -. 5195-00 f 

. 1847 - 00 ) 

. 1288+01 

( . 6104-01 f 

. 7509 - 00 ) 

. 2859+01 

( -. 5244-00 f 

. 1666 - 00 ) 

. 1319+01 

( . 5224-01 f 

. 7472 - 00 ) 

, 2890+01 

( -. 5288-00 f 

. 1484 - 00 ) 

, 1351+01 

( . 4265-01 f 

. 7431 - 00 ) 

. 2922+01 

( -. 5327-00 f 

. 1301 - 00 ) 

, 1382+01 

( . 3234-01 f 

. 7385 - 00 ) 

. 2953+01 

( -. 5361-00 f 

. 1117 + 00 ) 

, 1414+01 

( , 2137-01 f 

, 7336 - 00 ) 

, 2985+01 

( -. 5389-00 f 

. 9318 - 01 ) 

, 1445+01 

t . 9794-02 f 

. 7282 - 00 ) 

. 3016+01 

( -. 5413-00 f 

. 7463 - 01 ) 

, 1477+01 

( -. 2322-02 f 

. 7225 - 00 ) 

. 3047+01 

( -. 5431-00 f 

. 5602 - 01 ) 

, 1508+01 

( -. 1493-01 f 

. 7165 - 00 ) 

. 3079+01 

( -. 5444-00 f 

. 3737 - 01 ) 

, 1539+01 

( -. 2796-01 f 

, 7101 - 00 ) 

, 3110+01 

( -. 5452-00 f 

. 1869 - 01 ) 

, 1571+01 

( -. 4138-01 f 

, 7034 - 00 ) 

. 3142+01 

( -. 5455-00 f 

. 6097 - 15 ) 



TABLE AH. - LINE PRINTOUTS OF ADAMS -BASH FORTH PREDICTORS - Concluded 

(d) k = 4, d = 4 


ALPHA < 4) = 
ALPHA ( 3) = 
ALPHA < 2) = 
ALPHA ( 1) = 
ALPHA ( 0) = 


THETA 

,0000 

,3142-01 

.6283-01 

,9425-01 

.1257-00 

,1571-00 

,1885-00 

.2199-00 

,2513-00 

,2827-00 

.3142-00 

,3456-00 

,3770-00 

,4084-00 

.4398-00 

,4712-00 

.5027-00 

.5341-00 

,5655-00 

,5969-00 

,6283-00 

,6597-00 

.6912-00 

.7226-00 

,7540-00 

,7854-00 

,8168-00 

.8482-00 

.8796-00 

.9111-00 

.9425-00 

.9739-00 

.1005+01 

,1037+01 

.1068+01 

, 1100+01 

.1131+01 

.1162+01 

.1194+01 

,1225+01 

•1257+01 

,1288+01 

.1319+01 

,1351+01 

,1382+01 

.1414+01 

.1445+01 

.1477+01 

,1508+01 

.1539+01 

•1571+01 


ADAMS BASHFORTH, K=4* D=4 


• 1000000000000000+01 

-. 1000000000000000+01 

,0000000000000000 

.0000000000000000 

.QOOQOOOOOOOOOOOO 


BETA ( 4) 
BETAt 3) 
BETA ( 2) 
BETA ( 1) 
BETAt 0) 


.0000000000000000 
,2291666666666666+01 
-.2458333333333332+01 
. 1541666666666666+01 
-.3750000000000000+00 


HBAR THETA HBAR 



.0000 » 

-.0000 ) 

•1571+01 

( .2720-00 » 

,5779-00) 


.5205-09 , 

.3142-01) 

.1602+01 

( ,2524-00 * 

.5650-00) 


.3326-07 * 

,6283-01) 

,1634+01 

{ .2327-00 r 

.5529-00) 


.3779-06 #> 

.9425-01) 

.1665+01 

( .2130-00 > 

,5415-00) 


.2116-05 * 

,1257-00) 

.1696+01 

( ,1933-00 9 

.5306-00) 


,8038-05 p 

.1571-00) 

,1728+01 

( .1738-00 » 

,5202-00) 


.2388-04 # 

•1886-00) 

•1759+01 

( .1545-00 f 

.5102-00) 


•5985-04 f 

,2201-00) 

.1791+01 

( .1354-00 » 

.5005-00) 


.1325-03 , 

.2516-00) 

,1822+01 

( .1166+00 » 

.4911-00) 


.2666-03 p 

.2833-00) 

,1854+01 

( .9814-01 * 

.4818-00) 


.4978-03 p 

.3151-00) 

.1885+01 

( .7997-01 f 

.4726-00) 


.8747-03 p 

.3470-00) 

.1916+01 

t .6214-01 f 

.4635-00) 


.1462-02 * 

.3792-00) 

.1948+01 

( .4467-01 9 

.4544-00) 


,2341-02 r 

.4116-00) 

.1979+01 

( .2756-01 9 

.4453-00) 


.3617-02 r 

.4442-00) 

.2011+01 

{ .1083-01 9 

.4361-00) 


.5418-02 p 

.4771-00) 

,2042+01 

( -.5516-02 r 

.4268-00) 


.7898-02 » 

.5104-00) 

.2073+01 

( -.2147-01 » 

.4175-00) 


.1124-01 9 

.5439-00) 

,2105+01 

( -.3703-01 9 

.4080-00) 


.1565-01 9 

,5778-00) 

,2136+01 

( -.5218-01 » 

.3983-00) 


.2137-01 » 

.6118-00) 

,2168+01 

( -.6693-01 9 

.3885-00) 


.2867-01 f 

,6459-00) 

.2199+01 

( -.8126-01 » 

.3786-00) 


.3782-01 * 

.6799-00) 

,2231+01 

( -.9518-01 f 

.3684-00) 


.4910-01 9 

•7136-00) 

.2262+01 

( -.1087+00 » 

.3581-00) 


.6280-01 9 

.7465-00) 

.2293+01 

( -.1217+00 > 

.3475-00) 


.7914-01 » 

.7783-00) 

.2325+01 

( -.1344-00 9 

.3368-00) 


.9828-01 » 

.8084-00) 

.2356+01 

( -.1466-00 9 

.3259-00) 


.1203+00 * 

.8362-00) 

.2388+01 

( -.1584-00 9 

.3147-00) 


.1450-00 r 

,8609-00) 

.2419+01 

( -.1697-00 9 

.3034-00) 


.1722-00 r 

,8817-00) 

.2450+01 

( -.1806-00 9 

.2918-00) 


.2013-00 9 

.8981-00) 

,2482+01 

( -.1910-00 9 

.2801-00) 


.2317-00 r 

.9094-00) 

.2513+01 

t -.2010-00 r 

.2682-00) 


.2624-00 9 

.9151-00) 

.2545+01 

( -.2105-00 * 

.2560-00) 


.2924-00 r 

.9152-00) 

.2576+01 

( -.2195-00 » 

.2437-00) 


.3207-00 » 

.9096-00) 

,2608+01 

( -.2281-00 r 

.2312-00) 


.3462-00 f 

.8989-00) 

.2639+01 

t -.2363-00 > 

.2186-00) 


.3682-00 9 

.8836-00) 

.2670+01 

t -.2439-00 r 

.2058-00) 


.3859-00 » 

.8647-00) 

.2702+01 

( -.2511-00 * 

.1928-00) 


.3991-00 » 

.8429-00) 

.2733+01 

( -.2578-00 r 

.1797-00) 


.4076-00 9 

.6194-00) 

.2765+01 

( -.2640-00 # 

.1664-00) 


,4115-00 9 

.7948-00) 

•2796+01 

( -.2697-00 * 

.1530-00) 


.4112-00 f 

.7700-00) 

.2827+01 

( -.2749-00 9 

,1395-00) 


.4071-00 , 

.7456-00) 

.2859+01 

( -.2797-00 > 

.1258-00) 


.3996-00 9 

.7220-00) 

.2890+01 

( -.2839-00 9 

.1121+00) 


.3893-00 9 

.6994-00) 

,2922+01 

( -.2877-00 r 

.9831-01) 


•3766-00 9 

,6782-00) 

,2953+01 

( -.2910-00 9 

.8441-01) 


.3619-00 , 

.6582-00) 

,2985+01 

t -.2937-00 r 

.7045-01) 


,3457-00 , 

.6397-00) 

•3016+01 

( -.2960-00 r 

.5643-01) 


.3284-00 r 

,6225-00) 

.3047+01 

( -.2977-00 r 

.4236-01) 


•3101-00 , 

•6065-00) 

,3079+01 

( -.2990-00 9 

.2826-01) 


.2913-00 > 

•5917-00) 

.3110+01 

( -.2997-00 r 

.1414-01) 


.2720-00 , 

.5779-00) 

•3142+01 

( -.3000-00 9 

.4606-15) 


TABLE AIH. - LINE PRINTOUTS OF ADAMS -MOULTON CORRECTORS 


ALPHA < 1) = 
ALPHA t 0) = 


THETA 

.0000 

.3142-01 

,6283-01 

,9425-01 

.1257-00 

.1571-00 

.1885-00 

,2199-00 

,2513-00 

.2827-00 

,3142-00 

.3456-00 

.3770-00 

.4084-00 

.4398-00 

.4712-00 

.5027-00 

.5341-00 

.5655-00 

,5969-00 

.6283-00 

.6597-00 

.6912-00 

.7226-00 

.7540-00 

.7854-00 

,8168-00 

.8482-00 

,8796-00 

.9111-00 

.9425-00 

,9739-00 

.1005+01 

.1037+01 

.1068+01 

. 1100+01 

,1131+01 

.1162+01 

.1194+01 

.1225+01 

.1257+01 

.1286+01 

.1319+01 

.1351+01 

.1382+01 

.1414+01 

.1445+01 

.1477+01 

.1508+01 

,1539+01 

,1571+01 


(a) k = 1, d = 1 


ADAMS MOULTON K=1p D=1 


• 1000000000000000+01 

-.1000000000000000+01 


BETA< 1) = .1000000000000000+01 

BETA ( 0) = .0000000000000000 


HBAR THETA HBAR 



.0000 p 

-.0000 ) 

,1571+01 

( .1000+01 9 

,1000+01) 


.4934-03 p 

.3141-01) 

.1602+01 

( .1031+01 p 

.9995-00) 


,1973-02 p 

.6279-01) 

.1634+01 

< .1063+01 p 

,9980-00) 


.4438-02 f 

.9411-01) 

.1665+01 

( .1094+01 p 

.9956-00) 


.7885-02 p 

.1253-00) 

,1696+01 

( .1125+01 f 

.9921-00) 


.1231-01 t 

.1564-00) 

.1728+01 

( .1156+01 p 

.9877-00) 


.1771-01 p 

.1874-00) 

.1759+01 

( .1187+01 9 

.9823-00) 


.2408-01 r 

.2181-00) 

.1791+01 

< .1218+01 p 

.9759-00) 


.3142-01 f 

.2487-00) 

.1822+01 

( .1249+01 * 

.9686-00) 


,3971-01 r 

,2790-00) 

.1854+01 

( .1279+01 * 

.9603-00) 


,4894-01 p 

.3090-00) 

.1885+01 

( .1309+01 p 

.9511-00) 


.5912-01 p 

.3387-00) 

.1916+01 

( .1339+01 p 

.9409-00) 


.7022-01 * 

.3601-00) 

.1948+01 

< ,1368+01 p 

.9298-00) 


.8225-01 * 

,3971-00) 

.1979+01 

< ,1397+01 p 

.9178-00) 


.9517-01 p 

.4258-00) 

,2011+01 

( .1426+01 p 

.9048-00) 


.1090+00 p 

.4540-00) 

.2042+01 

( .1454+01 p 

.8910-00) 


.1237+00 p 

.4818-00) 

,2073+01 

( .1482+01 p 

.8763-00) 


,1393-00 p 

.5090-00) 

.2105+01 

( .1509+01 p 

.8607-00) 


.1557-00 t 

.5358-00) 

,2136+01 

( .1536+01 p 

.8443-00) 


.1729-00 r 

,5621-00) 

.2168+01 

( ,1562+01 p 

.8271-00) 


.1910-00 p 

,5878-00) 

,2199+01 

( .1588+01 p 

.8090-00) 


.2098-00 9 

,6129-00) 

.2231+01 

( .1613+01 p 

.7902-00) 


.2295-00 * 

.6374-00) 

.2262+01 

( .1637+01 p 

.7705-00) 


,2499-00 p 

.6613-00) 

,2293+01 

( ,1661+01 p 

.7501-00) 


.2710-00 9 

.6845-00) 

.2325+01 

( .1685+01 p 

.7290-00) 


.2929-00 9 

.7071-00) 

.2356+01 

( .1707+01 p 

.7071-00) 


,3155-00 9 

.7290-00) 

.2388+01 

( .1729+01 p 

,6845-00) 


,3387-00 * 

.7501-00) 

.2419+01 

( .1750+01 p 

.6613-00) 


.3626-00 p 

.7705-00) 

.2450+01 

( .1771+01 p 

.6374-00) 


,3871-00 p 

.7902-00) 

,2482+01 

( .1790+01 p 

.6129-00) 


.4122-00 9 

,8090-00) 

,2513+01 

( .1809+01 p 

.5878-00) 


.4379-00 p 

.8271-00) 

,2545+01 

( .1827+01 p 

.5621-00) 


,4642-00 * 

,8443-00) 

,2576+01 

( .1844+01 p 

.5358-00) 


.4910-00 p 

,8607-00) 

,2608+01 

( ,1861+01 p 

.5090-00) 


,5182-00 p 

.8763-00) 

.2639+01 

( .1876+01 p 

.4818-00) 


,5460-00 9 

.8910-00) 

.2670+01 

( .1891+01 p 

.4540-00) 


,5742-00 p 

.9048-00) 

.2702+01 

( .1905+01 p 

.4258-00) 


,6029-00 p 

.9178-00) 

.2733+01 

( .1918+01 p 

.3971-00) 


.6319-00 p 

.9298-00) 

.2765+01 

{ .1930+01 p 

,3681-00) 


.6613-00 9 

.9409-00) 

.2796+01 

( .1941+01 p 

.3387-00) 


.6910-00 * 

,9511-00) 

.2827+01 

( .1951+01 p 

,3090-00) 


.7210-00 p 

.9603-00) 

,2859+01 

( .1960+01 p 

.2790-00) 


,7513-00 9 

•9606-00) 

.2890+01 

( .1969+01 p 

.2487-00) 


.7819-00 p 

,9759-00) 

.2922+01 

( .1976+01 p 

.2181-00) 


.8126-00 9 

.9823-00) 

,2953+01 

( .1982+01 p 

,1874-00) 


.8436-00 * 

.9877-00) 

.2985+01 

( .1988+01 p 

.1564-00) 


,8747-00 p 

,9921-00) 

.3016+01 

( .1992+01 p 

.1253-00) 


.9059-00 p 

.9956-00) 

.3047+01 

( .1996+01 p 

.9411-01) 


,9372-00 9 

.9980-00) 

.3079+01 

{ .1998+01 p 

.6279-01) 


.9686-00 p 

.9995-00) 

.3110+01 

( .2000+01 p 

.3141-01) 


.1000+01 9 

.1000+01) 

.3142+01 

( .2000+01 p 

.9992-15) 



TABLE AID. - LINE PRINTOUTS OF ADAMS -MOULTON CORRECTORS - Continued 

(b) k = 1, d = 2 


ALPHA ( 1) = 
ALPHA ( 0) = 


THETA 

.0000 

.3142-01 

.6283-01 

.9425-01 

.1257-00 

.1571-00 

.1885-00 

.2199-00 

.2513-00 

.2827-00 

.3142-00 

,3456-00 

.3770-00 

.4084-00 

.4398-00 

,4712-00 

.5027-00 

.5341-00 

.5655-00 

.5969-00 

.6283-00 

.6597-00 

.6912-00 

.7226-00 

.7540-00 

.7854-00 

.8168-00 

,8482-00 

,8796-00 

,9111-00 

.9425-00 

.9739-00 

.1005+01 

.1037+01 

.1068+01 

. 1100+01 

.1131+01 

.1162+01 

.1194+01 

.1225+01 

,1257+01 

,1288+01 

.1319+01 

,1351+01 

.1382+01 

.1414+01 

.1445+01 

.1477+01 

.1508+01 

.1539+01 


ADAMS MOULTON K=1 f 0=2 


,1000000000000000+01 BETAt 1) = .5000000000000000+00 

-.1000000000000000+01 BETA ( 0) = .5000000000000000+00 


HBAR THETA HBAR 



•0000 r 

-.0000 ) 

.1571+01 

{ -.2220-15 9 

•2000+01) 


-.6984-16 f 

,3142-01) 

•1602+01 

{ -.2866-15 r 

•2064+01) 


-.1189-15 f 

.6285-01) 

.1634+01 

( -.2962-15 f 

.2130+01) 


-.2260-16 f 

.9432-01) 

.1665+01 

( -.4902-15 f 

.2198+01) 


-.7924-16 r 

.1258-00) 

.1696+01 

{ -.6347-16 f 

,2269+01) 


.2618-16 f 

.1574-00) 

,1728+01 

( -.6581-16 r 

.2342+01) 


,2450-16 * 

.1891-00) 

.1759+01 

( .0000 9 

.2418+01) 


.8779-17 f 

.2208-00) 

.1791+01 

( -.7100-16 9 

.2496+01) 


-.7050-16 t 

.2527-00) 

.1822+01 

( -.2217-15 f 

,2578+01) 


-.6371-16 f 

.2846-00) 

.1854+01 

( -.3850-15 9 

.2664+01) 


-.8891-16 r 

.3168-00) 

,1885+01 

( -.4820-15 9 

.2753+01) 


-.8223-16 * 

.3491-00) 

.1916+01 

( -.4197-15 f 

.2846+01) 


-.1223-15 t 

.3815-00) 

.1948+01 

( -.5271-15 » 

.2943+01) 


-.1013-15 * 

.4142-00) 

.1979+01 

( -.3683-15 f 

.3045+01) 


-.2914-16 * 

.4471-00) 

.2011+01 

( .0000 9 

.3151+01) 


-.8807-16 r 

.4802-00) 

.2042+01 

( -.1017-15 f 

.3264+01) 


.7396-17 f 

.5135-00) 

.2073+01 

< -.4285-15 9 

,3382+01) 


-.1492-16 r 

.5471-00) 

.2105+01 

( -.3392-15 * 

.3506+01) 


-.6020-16 f 

.5811-00) 

.2136+01 

( -.4784-15 9 

•3638+01) 


-.1671-15 » 

.6153-00) 

.2168+01 

( -.1014-14 f 

.3777+01) 


-.1381-15 t 

,6498-00) 

.2199+01 

( -.6733-15 9 

.3925+01) 


-.1240-15 f 

.6848-00) 

.2231+01 

( -.7170-15 f 

.4083+01) 


-.1254-15 9 

.7200-00) 

.2262+01 

( -.6124-15 f 

.4250+01) 


-.6344-16 9 

.7557-00) 

,2293+01 

( -.1639-15 , 

.4430+01) 


-.1605-15 f 

.7919-00) 

,2325+01 

( -.1056-14 9 

.4622+01) 


-.1301-15 9 

.8284-00) 

.2356+01 

{ -.9476-15 » 

.4828+01) 


-.1977-15 f 

.8655-00) 

.2388+01 

( -.8193-15 » 

.5051+01) 


-.1671-15 f 

.9030-00) 

.2419+01 

( -.6664-15 * 

.5293+01) 


-.6780-16 f 

,9411-00) 

.2450+01 

( -.3628-15 9 

.5555+01) 


-.1377-15 f 

.9798-00) 

•2482+01 

( -.5291-15 , 

.5842+01) 


-.2447-15 9 

.1019+01) 

,2513+01 

( -.1163-14 r 

.6155+01) 


-.1777-15 , 

,1059+01) 

,2545+01 

( -.1605-14 > 

.6501+01) 


-.1807-15 9 

,1100+01) 

.2576+01 

( -.1605-14 9 

.6884+01) 


-.7357-16 * 

.1141+01) 

.2608+01 

( -.5979-15 > 

.7311+01) 


-.1873-15 * 

.1183+01) 

.2639+01 

( -.1234-14 * 

.7789+01) 


-.3818-16 f 

,1226+01) 

.2670+01 

{ .1273-15 r 

.8331+01) 


-.3893-16 f 

.1269+01) 

.2702+01 

( -.7291-15 , 

.8947+01) 


-.1192-15 f 

.1314+01) 

.2733+01 

( -.2531-14 * 

•9658+01) 


-.1623-15 f 

.1359+01) 

.2765+01 

( -.3755-14 » 

.1048+02) 


-.1659-15 f 

.1406+01) 

.2796+01 

( -.3169-14 r 

.1146+02) 


-.4241-15 t 

.1453+01) 

.2827+01 

( -.3403-14 t 

.1263+02) 


-.2170-15 $ 

.1502+01) 

.2859+01 

( -.4544-14 f 

.1405+02) 


-.1334-15 9 

,1551+01) 

.2890+01 

( -.3313-14 f 

.1583+02) 


-.4557-16 t 

.1602+01) 

.2922+01 

( -.6051-14 f 

.1812+02) 


.0000 * 

.1655+01) 

.2953+01 

( -.5876-14 f 

.2116+02) 


-.4800-16 f 

.1708+01) 

.2985+01 

< -.4932-14 f 

.2541+02) 


-.9866-16 * 

.1763+01) 

.3016+01 

( -.1452-13 f 

.3179+02) 


-.4059-15 f 

.1820+01) 

.3047+01 

( -.2013-13 f 

.4241+02) 


-.3656-15 9 

.1878+01) 

.3079+01 

( -.6066-13 f 

.6364+02) 


-.2691-15 f 

.1938+01) 

.3110+01 

( -.5632-12 f 

.1273+03) 



TABLE AID. - LINE PRINTOUTS OF ADAMS -MOULTON CORRECTORS - Continued 


ALPHA ( 2) = 
ALPHA ( 1) = 
ALPHA { 0) = 


THETA 

,0000 
,3142-01 
,6283-01 
,9425-01 
,1257-00 
,1571-00 
,1885-00 
,2199-00 
.2513-00 
.2827-00 
.3142-00 
.3456-00 
.3770-00 
,4084-00 
.4398-00 
.4712-00 
,5027-00 
.5341-00 
.5655-00 
.5969-00 
,6283-00 
.6597-00 
.6912-00 
,7226-00 
.7540-00 
,7854-00 
,8168-00 
.8482-00 
.8796-00 
,9111-00 
.9425-00 
,9739-00 
.1005+01 
,1037+01 
,1068+01 
. 1100+01 
.1131+01 
,1162+01 
.1194+01 
, 1$25+01 
.1257+01 
.1288+01 
.1319+01 
.1351+01 
.1382+01 
.1414+01 
,1445+01 
.1477+01 
,1508+01 
.1539+01 
,1571+01 


(c) k = 2, d = 3 


ADAMS MOULTON K=2f D=3 


. 1000000000000000+01 

-. 1000000000000000+01 

,0000000000000000 


BETAt 2) = ,4166666666666666+00 
BETA ( 1) = .6666666666666665+00 
BETA( 0) = -.8333333333333331-01 


HBAR THETA HBAR 



.0000 r 

-.0000 ) 

.1571+01 

( -.2400-00 f 

.1680+01) 


-.4058-07 f 

.3142-01) 

.1602+01 

{ -.2606-00 f 

,1722+01) 


-.6492-06 r 

.6283-01) 

.1634+01 

( -.2827-00 f 

.1764+01) 


-.3285-05 f 

,9425-01) 

,1665+01 

{ -.3063-00 f 

.1807+01) 


-.1038-04 r 

.1257-00) 

.1696+01 

( -.3316-00 f 

.1851+01) 


-.2532-04 f 

,1571-00) 

.1728+01 

( -.3586-00 f 

.1895+01) 


-.5245-04 t 

,1885-00) 

.1759+01 

( -.3876-00 f 

.1941+01) 


-.9706-04 f 

.2199-00) 

.1791+01 

< -.4185-00 f 

.1987+01) 


-.1654-03 f 

.2513-00) 

,1822+01 

< -.4517-00 f 

.2034+01) 


-.2646-03 f 

.2828-00) 

.1854+01 

( -.4872-00 f 

.2082+01) 


-.4027-03 r 

.3142-00) 

,1885+01 

( -.5251-00 f 

.2130+01) 


-.5886-03 f 

,3456-00) 

,1916+01 

( -.5658-00 f 

.2180+01) 


-.8322-03 » 

.3771-00) 

.1948+01 

( -.6093-00 f 

.2230+01) 


-.1144-02 f 

.4086-00) 

.1979+01 

( -.6559-00 f 

.2281+01) 


-.1536-02 f 

.4401-00) 

,2011+01 

( -.7058-00 f 

.2333+01) 


-.2020-02 f 

,4716-00) 

.2042+01 

( -.7593-00 f 

.2385+01) 


-.2609-02 f 

,5031-00) 

,2073+01 

( -.8166-00 f 

,2439+01) 


-.3318-02 f 

,5347-00) 

.2105+01 

< -.8781-00 f 

.2492+01) 


-.4161-02 f 

,5663-00) 

.2136+01 

( -.9441-00 f 

.2547+01) 


-.5153-02 f 

,5980-00) 

,2168+01 

( -.1015+01 f 

.2601+01) 


-.6312-02 f 

,6297-00) 

.2199+01 

( -.1091+01 f 

,2656+01) 


-.7654-02 f 

,6615-00) 

.2231+01 

( -.1172+01 f 

.2711+01) 


-.9196-02 f 

.6934-00) 

.2262+01 

( -.1260+01 f 

.2765+01) 


-.1096-01 f 

.7254-00) 

.2293+01 

( -.1354+01 f 

.2819+01) 


-.1296-01 f 

.7574-00) 

.2325+01 

( -.1455+01 f 

.2872+01) 


-.1522-01 f 

.7896-00) 

.2356+01 

{ -.1563+01 f 

.2923+01) 


-.1776-01 f 

.8219-00) 

,2388+01 

{ -.1679+01 f 

.2973+01) 


-.2061-01 f 

.8543-00) 

.2419+01 

< -.1804+01 f 

.3019+01) 


-.2378-01 f 

.8869-00) 

,2450+01 

< -.1938+01 f 

.3063+01) 


-.2730-01 f 

.9196-00) 

.2482+01 

( -.2082+01 f 

.3102+01) 


-.3120-01 f 

.9525-00) 

.2513+01 

{ -.2236+01 f 

,3136+01) 


-.3550-01 f 

,9856-00) 

,2545+01 

( -.2400+01 f 

.3163+01) 


-.4023-01 f 

,1019+01) 

.2576+01 

( -.2575+01 f 

.3182+01) 


-.4541-01 f 

.1053+01) 

,2608+01 

< -.2761+01 f 

.3191+01) 


-.5108-01 f 

,1086+01) 

.2639+01 

( -.2959+01 f 

.3189+01) 


-.5727-01 f 

,1120+01) 

.2670+01 

< -.3167+01 f 

.3173+01) 


-.6401-01 f 

.1155+01) 

.2702+01 

( -.3386+01 f 

.3141+01) 


-.7134-01 f 

.1190+01) 

.2733+01 

( -.3615+01 f 

.3091+01) 


-.7929-01 f 

.1225+01) 

.2765+01 

( -.3852+01 f 

.3019+01) 


-.8790-01 f 

.1260+01) 

.2796+01 

{ -.4095+01 f 

,2924+01) 


-.9722-01 f 

,1296+01) 

.2827+01 

( -.4342+01 f 

.2802+01) 


-.1073+00 f 

.1332+01) 

.2859+01 

( -.4589+01 f 

,2652+01) 


-.1181+00 f 

.1369+01) 

.2890+01 

( -.4832+01 f 

,2471+01) 


-.1298-00 f 

.1406+01) 

.2922+01 

( -.5067+01 f 

,2258+01) 


-.1424-00 f 

.1443+01) 

.2953+01 

( -.5288+01 f 

.2013+01) 


-.1559-00 f 

,1481+01) 

.2985+01 

( -.5488+01 f 

.1736+01) 


-.1705-00 » 

,1520+01) 

.3016+01 

( -.5663+01 f 

.1429+01) 


-.1861-00 f 

.1559+01) 

.3047+01 

( -.5806+01 f 

.1097+01) 


-.2028-00 f 

.1599+01) 

,3079+01 

( -.5912+01 f 

.7438-00) 


-.2208-00 f 

.1639+01) 

.3110+01 

( -.5978+01 f 

.3757-00) 


-.2400-00 f 

,1680+01) 

,3142+01 

( -.6000+01 f 

.1199-13) 



TABLE Am. - LINE PRINTOUTS OF ADAMS -MOULTON CORRECTORS - Concluded 

(d) k = 3, d = 4 


ADAMS MOULTON K=3» D=4 


ALPHA { 3) = ,1000000000000000+01 

ALPHA ( 2) = -.1000000000000000+01 
ALPHA! 1) = .0000000000000000 

ALPHA ( 0) = ,0000000000000000 


BETA ( 3) = .3750000000000000+00 
BETAt 2) = .7916666666666665+00 
BETA ( 1) = -.2083333333333332+00 
BETA ( 0) = .4166666666666666-01 


THETA 

HBAR 

THETA 

HBAR 


• 0000 

( ,0000 * 

-.0000 ) 

.1571+01 

( -.1846-00 

t 

.1477+01) 

.3142-01 

( -.2003-10 » 

.3142-01) 

.1602+01 

( -.2037-00 

r 

.1502+01) 

.6283-01 

( -.1281-08 • 

.6283-01) 

,1634+01 

( -.2242-00 

f 

.1527+01) 

,9425-01 

( -.1458-07 9 

,9425-01) 

.1665+01 

( -.2462-00 

t 

.1552+01) 

.1257-00 

( -.8182-07 t 

.1257-00) 

,1696+01 

{ -.2697-00 

t 

.1576+01) 

.1571-00 

C -.3117-06 • 

.1571-00) 

.1728+01 

( -.2949-00 

t 

.1600+01) 

,1885-00 

( -.9289-06 » 

.1885-00) 

.1759+01 

< -.3218-00 

t 

.1623+01) 

.2199-00 

( -.2337-05 r 

.2199-00) 

.1791+01 

( -.3505-00 

f 

.1647+01) 

.2513-00 

( -.5194-05 * 

.2513-00) 

,1822+01 

( -.3810-00 

9 

.1669+01) 

.2827-00 

< -.1050-04 * 

,2827-00) 

,1854+01 

< -.4135-00 

9 

.1692+01) 

.3142-00 

( -.1969-04 r 

.3141-00) 

,1885+01 

< -.4479-00 

9 

.1714+01) 

.3456-00 

( -.3476-04 » 

.3455-00) 

,1916+01 

( -.4844-00 

9 

.1735+01) 

.3770-00 

< -.5835-04 » 

.3768-00) 

,1948+01 

( -.5231-00 

9 

.1755+01) 

.4084-00 

( -.9389-04 t 

.4081-00) 

.1979+01 

{ -.5641-00 

9 

.1775+01) 

.4398-00 

( -.1457-03 , 

.4394-00) 

.2011+01 

( -.6074-00 

9 

.1793+01) 

,4712-00 

( -.2193-03 » 

.4707-00) 

.2042+01 

( -.6530-00 

9 

.1811+01) 

.5027-00 

( -.3212-03 , 

.5019-00) 

.2073+01 

( -.7012-00 

9 

.1828+01) 

,5341-00 

( -.4592-03 * 

.5330-00) 

.2105+01 

( -.7519-00 

9 

.1843+01) 

.5655-00 

( -.6427-03 r 

.5641-00) 

,2136+01 

( -.8052-00 

9 

.1857+01) 

.5969-00 

< -.8826-03 t 

.5952-00) 

,2168+01 

( -.8612-00 

9 

.1869+01) 

.6283-00 

( -.1192-02 > 

.6261-00) 

.2199+01 

( -.9199-00 

9 

.1880+01) 

.6597-00 

( -.1584-02 * 

.6570-00) 

.2231+01 

( -.9815-00 

9 

.1888+01 ) 

.6912-00 

( -.2076-02 r 

.6877-00) 

.2262+01 

( -.1046+01 

9 

.1894+01) 

.7226-00 

( -.2685-02 » 

.7183-00) 

.2293+01 

( -.1113+01 

9 

.1897+01) 

.7540-00 

( -.3433-02 » 

.7488-00) 

.2325+01 

{ -.1183+01 

9 

.1898+01) 

,7854-00 

( -.4341-02 * 

.7792-00) 

.2356+01 

( -.1256+01 

9 

.1896+01) 

.8168-00 

( -.5434-02 # 

.8095-00) 

.2388+01 

( -.1331+01 

9 

.1890+01) 

.8482-00 

( -.6738-02 r 

.8395-00) 

.2419+01 

( -.1410+01 

9 

,1880+01) 

.8796-00 

( -.8282-02 » 

.8695-00) 

,2450+01 

( -.1490+01 

9 

.1867+01) 

,9111-00 

( -.1010-01 » 

,8992-00 ) 

,2482+01 

( -.1573+01 

9 

.1848+01) 

,9425-00 

( -.1222-01 » 

,9288-00) 

,2513+01 

( -.1659+01 

9 

.1825+01) 

,9739-00 

t -.1467-01 , 

.9581-00) 

.2545+01 

( -.1746+01 

9 

.1797+01) 

.1005+01 

( -.1751-01 » 

.9873-00) 

.2576+01 

( -.1835+01 

9 

. 1762+01 ) 

.1037+01 

( -.2076-01 . 

.1016+01) 

.2608+01 

{ -.1925+01 

9 

.1722+01) 

,1068+01 

! -.2446-01 > 

,1045+01) 

.2639+01 

( -.2016+01 

9 

.1676+01) 

,1100+01 

( -.2866-01 > 

,1074+01) 

.2670+01 

( -.2107+01 

9 

.1622+01 ) 

,1131+01 

( -.3340-01 » 

,1102+01) 

,2702+01 

( -.2198+01 

9 

.1562+01) 

,1162+01 

( -.3873-01 » 

.1130+01) 

,2733+01 

{ -.2288+01 

9 

.1494+01) 

.1194+01 

t -.4469-01 # 

,1158+01) 

.2765+01 

( -.2376+01 

9 

.1418+01) 

,1225+01 

( -.5134-01 » 

.1186+01) 

,2796+01 

( -.2461+01 

9 

. 1335+01 ) 

,1257+01 

1 -.5872-01 r 

,1213+01) 

.2827+01 

( -.2544+01 

9 

.1244+01) 

,1288+01 

( -.6688-01 » 

.1241+01) 

,2859+01 

( -.2622+01 

9 

.1146+01) 

.1319+01 

( -.7588-01 f 

.1268+01) 

,2890+01 

< -.2695+01 

9 

,1040+01) 

.1351+01 

< -.8576-01 t 

.1295+01) 

.2922+01 

( -.2762+01 

9 

.9278-00) 

,1382+01 

( -.9659-01 » 

.1321+01) 

,2953+01 

( -.2822+01 

9 

.8087-00) 

.1414+01 

( -.1084+00 . 

.1348+01) 

,2985+01 

( -.2875+01 

9 

,6836-00) 

,1445+01 

t -.1213+00 . 

,1374+01) 

.3016+01 

( -.2919+01 

9 

,5535-00) 

,1477+01 

( -.1353-00 » 

.1400+01) 

.3047+01 

( -.2954+01 

9 

.4190-00) 

,1508+01 

( -.1505-00 $ 

.1426+01) 

.3079+01 

( -.2979+01 

9 

,2812-00) 

.1539+01 

< -.1669-00 , 

,1452+01) 

,3110+01 

( -.2995+01 

9 

.1412-00) 

,1571+01 

( -.1846-00 r 

.1477+01) 

.3142+01 

( -.3000+01 

9 

.4673-14) 




-3 - 2-10123 

(a) k = 1, d = 1 

Figure Al. - FORTRAN IV subroutine STBL1 plots of Adams-Bashforth predictors. 
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-3 -2-10123 

(d) k = 4, d = 4 
Figure Al. - Concluded. 
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(a) k = 1, d = 1 

Figure A2. - FORTRAN IV subroutine STBLl plots of Adams-Moulton correctors. 
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APPENDIX B 


FORTRAN IV SUBROUTINE STBL2 FOR PREDICTOR- CORRECTOR 

MULTISTEP METHODS 

The subroutine STBL2 computes the boundary of the region of stability for 
predictor-corrector multistep methods for which the corrector is applied only once. 

In order to use the subroutine STBL2, the following calling sequence is necessary: 

CALL STBL2 (AA1, BB1, AA2, BB2, K, N, TITLE) 

Where the predictor-corrector multistep method is 

Predictor: « k y n+k + « k _ 1 y n4k _ 1 + ••• + <V n = h (\-l y 'n+k-l + ”* + Vn) 

Corrector: y R y n+k + 7^7^^ + • • • + r 0 y n = h ( 6 k y ’n+k + 6 k-l y 'n+k-l 

+ . . . + 6 ^’^ 

and 

AA1 is a double precision array with 25 locations such that AA1(1) = a k> 

AA1(2) = « k • • • , AAl(k + 1) = Oq. 

BBl is a double precision array with 25 locations such that BBl(l) = 0, 
BBl(2)-j9^ r , BBl(k + 1) = j8 0 . 

AA2 is a double precision array with 25 locations such that AA2(1) =y k , 

AA2(2) = r k _p • • • , AA2(k + 1) = y Q . 

BB2 is a double precision array with 25 locations such that BB2(1) = 6 k , 

BB2(2) = . . . , BB2(k + 1) = 6 Q . 

K is an integer giving the order of the method (use the largest order of the 
predictor or of the corrector). 

N is an integer giving the number of divisions of the interval (0, n) to be 
used in computing the boundary of the region of stability; 90 divisions 
are usually sufficient. 
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TITLE is an array of 5 locations containing BCD information to be used as a title 
for the printed and plotted output. TITLE may be passed to the subroutine 
as 30H (any 30 BCD characters). 

The subroutine CMPLX is used to perform the double precision complex arith- 
metic. 

STBL2 will compute the boundary of the region of stability and print the results 
on the line printer if 10 = 6 or on the 4020 output if 10 = 17; it will also plot the results 
on the 4020 plotter. 

Table BI is a listing of the FORTRAN IV subroutine STBL2. Subroutine STBL2 
has been used to calculate the region of stability for the Bashforth-Moulton predictor- 
corrector pairs of degree 1 to 4. Table BII gives the coefficients of Bashforth-Moulton 
methods and the numerical results of STBL2, while figure Bl displays the results 
graphically. The coefficients for the methods are available from NASA Manned 
Spacecraft Center, Houston, Texas. 
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TABLE BI. - LISTING OF FORTRAN IV SUBROUTINE STBL 2 FOR 
PREDICTOR -CORRECTOR MULTISTEP METHODS 


DIMENSION AA125) # Bb < 25 ) r aAX ( 2R ) f bBl ( 2S) » AA2 < 25 ) » bR2 ( 25) 
UOUBLE PKELlblON AA * bR r AA1 * BBl » Bb2 » AA? 

AA1 ( 1 ) = l .ObO 

AA1 (2) = -i.UDO 
A A2 ( 1 ) = 1.0U0 
AA2 ( 2 > = -1.0DU 
JjO 1U01 w=3*25 
AA1 ( J ) = O.ObO 
1001 A A 2 ( J ) =_ 0*0U0 
BR( 1 ) = b.uOU 
bB<2) = i.uDu 
CALL M0Vt(oB»RBl) 
bR(l) = 1.UD0 
bBl2) = U.UDQ 

CALL M0VL(bB*PB2> 

KS1EP = i 

CALL STBl 2< AA1 *Bbl»AA2»Bb2»KSTLP»P0» 

1 3uHbAbHhOkTH-MOULTON K = l* 0=1 ) 

Bb(i) = U.UDO 

Bb ( 2 ) = u.l5UOoOCUuOOOuOOObl 

BBJ3L =. -O.bOgOOOuOUOOODOaOOO 

CALL MOVE <DB, BBl) 

Bb ( 1) = 0-5u0 
PB(2) = 0-5un 
bB<3) = u.uDO 
CALL MOVE ( bB t BB2 ) 

ASTEH = 2 

CALL STBL21AA1 ,Bbl » AM2»Bb2»KSTLP*P0» 

1 30HBASHEOKTH-WOULTON K=2t 0=2 ) 

BB<1) = U.UDO 

Pb ( c ) = U . l9l6b6o6b666 b 6fo7Ul 
Rb<3> = -U.1333333333333333U1 

Bb(4) = li • 4 166bbbbo6b6 b 6t)7U0 

CALL MOVt(bBrPbl) 

Bb(l) = 0 • 4l6b666b6o6b6 b 67D-nO 
Pb(2) = 0*6 b 6b6b6b6b6b6b67D-nU 
BB ( 3 ) =-0.83333333333333330-01 
bB ( 4 ) = O.UDU 

CALL N>0VL(bB»Rb2) 

KS7EH = 3 

CALL STBl21Am 1 rBb.l r A /A ?»Bb2»KST£PrP0» 

1 3uHbAbHF- OkTH-N'OULT ON K = 3r 0=3 ) 

tiRll) = U.uDu 

Bb(^) = U.2291b6b6fo66666b7ui 

Bb(3) = -U. 2458333333333333b 1 . __ 

Bb { 4 ) = U. 1541b6b6b66666b7Dl 
Pb ( 5 ) = -u.o7b000U000U0(i0U0un 
CALL MOV/t ( bB t Bbl ) 

Bb ( 1 ) = 0.37500 

Bb ( 2 ) = U.791666b6b6b666b7U-on 

Pb(3> = 'U.2083333333333333U-00 

PtJ(4) = U • 4 Ib6b6bb6668b6b7u-U 1 
tiRl5> - U.UDU 
CALL M0Vb(bB.BB2) 
fVSTEP = 4 

CALL 5TBl 2 l AA1 » Bbl * A,\2 # Bb2 * KSTLP »9l)» 

I 3 UHb AiiHfc 0 K TH-jMO UL TON K=4» . 0=4 ) 

END 



TABLE BI. - LISTING OF FORTRAN IV SUBROUTINE STBL2 FOR 
PREDICTOR -CORRECTOR MULTISTEP METHODS - Continued 


SUBROUTINE MoVt(LjBrHbl) 
JImEnSIOi* oB ( ?S > t Bbl ( 2b ) 
uOUHLE PRECISION Bb'bPl 
uO 1U10 i=i*2S 
1UI0 uFUU)=BbU) 

RETURN 

lnu 


SUBROUTINE COMPLX (Ar IA, b» Tbt Kr Cr IC> 

INTEGER K 

DOUBLE PRECISION Ar IAr B# Ibr Cr ICr TEMP r DSQRT 

DOUBLE PRECISION TEMpC 

00 TO ( 1 f ** M » K 
C 1 = MPYr 2 = uvur .5 = SOT 

1 TFMPC = m ♦ b - IA * Tb 
iC = A * It> + IA * B 

C = TEMPO 
00 TO 10U 
21 1C = 0.0U0 
C - l.ODu 
00 TO 100 

2 CONTINUE 

TEMP = DABS ( b**2 + Ib**2> 

IF HEMP *LT, 1.0D-2u) GO TO 2 1 
1EMPC = IA * p + IA * IB) / TFMP 
1C = (1A * B - A * Itj) / TEMP 
C - IEMPC 
GO TO lOu 

b TEMP = DbQKT(DABS< A**? + IA**?)) 

TEMPO = uSoRT(DABS(T L MP + A) / 2.nD07 — 

MINUS = 1 

if (ia .lT. u.odo) Minus = -1 
ic = dsqrt i D abs ( temp - A) / 2.0D0) 
c - TEMPO 
1C = MlNUS*IC 

IF ((A .LT* U.ODO) *AuO* (IA .GT. 0*000) ) GO TO 90 
GO TU lOu 
V9 CONTINUL 
C = -C 
1C = -IC 
lou keiukn 

ENU 



TABLE BL - LISTING OF FORTRAN IV SUBROUTINE STBL2 FOR 
PREDICTOR -COLLECTOR MULTISTEP METHODS - Continued 


SUBROUTINE SI BL2( AAl *BB1 # AA2 r RB2 » K # Nr TITLE ) . 

DIMENSION AAl (25) rbBl(25) #AA2(25> »BB2 (25) r Hi (500) rh2(500) r 
I H315U0) fh4(B00) tTITLE(b) 

DIMENSION f H ( 500 > 

DIMENSION BCDX(12)» uCDY(12> 

DATA BCDa/12*FH / 

DA 1 A RCDI/A2 + 6H / 

DOUBLE PRECISION A A 1 » Rbl » A A2 » Rb2 » AK r A I »RR#bI*CR»Cl r Dk r D I r DC Os * 
i 0 S I J j » ACk » AC 1 1 BRk r bRI 9 DSr »DSI » HR » HI »_DL 

dourle precision Tutu* del 
10 = 6 
10 = 17 
K1=K+1 
1C=0 
NC-0 

UEL=3.14l5^2b535ti979^,nu/N 

1HLTA=0.U0U 

1 U20 IF (NC • bT • n) GO l 0 1030 
THtTA = NC*Dtl 
AR^O.OUO 
aJ^O.ODO 
bR=0. 0D0 
dT=O.ODO 
CR-0 • 0D0 
C I -0 • 0D0 
DO 1010 1 = 

JL=K-i *1 

DL=L ‘ " " " " 

IF (1 .Eu. 1) GO To m21 
aR— A h+bBl ( 1 ) *DCOS (dL*jHE1 A ) 

A I-A 1+BBi ( 1 >*OSIn(DL*tHFTA) 

bR = RR+(Rc52(I)-fcjB2(i)»AAl(I))*DLOS(Dt»THFTA) 

dl = BI+lBb2< I)-bH2<l)*AAl il) )*DSlN(b| *THFTA) 

1021 CONTINUE 

cr— ck-aa^ ( i ) *ocos ( UL« tHE7 A ) 

C I— C I —A A 2 ( 1 ) *DblN ( DL* tHE T A ) 

101U CONTINUE 

AR~AR*BB2 ( 1 ) 

Al-Al*bB«£f 1) 

CALL COVHLX ( 4*AR»4*Al *CRrC 1 r i r ACrt » ACI ) 

CALL CGMHLa ( bR r B 1 * bR * Rl * 1 »t>BRrbBl ) 

UR— BBR—ACR 
DI=BbI-AcI 

CALL C0MHLA(uR*01#l*R»nI*3»DSR»DSl ) 

UR=-bR+DSR 

uI = -oI+0Sl 

CALL COMHLA(DPrDl *2*A»»2^A1 f 2»HR»HI ) 

IC = IC + l 
HJ <IC> = HK 
H2UC) = HI 
IH <IC) = TNEIA 
UR = -BR-DSR 
UI = -til-Dbl 

CALL COMhLX < UR r Di r 2*ARt 2* Air 2, HRr HI ) 

IF (IC .LQ- 1) Go TO H4u 

D IF 1 = ( hi ( I c**l ) -HK ) **2 + ( H2 ( lr-1 ) —H I ) ** ? 

DIF 3 = (H3<lC-:l)-HK)**2+<H4(IC-n-Hn**2 

uiF2=(Hl(lC-l)-Hl( IC) )**2+lH2( IC-1 ) -h? ( IC))**? 

DIF4=(H31 IC-D-Hl ( iC))**2+( H4 ( IC — 1 ) **H2 ( I C ) >**? 

IF (IDlFo *LT ♦ DiFl ) .OR, (0IF2 .LT. nlF4))60 TO U4o 



TABLE BI. - LISTING OF FORTRAN IV SUBROUTINE STBL2 FOR 
PREDICTOR-COLLECTOR MULTISTEP METHODS - Concluded 


H34IU = HlllC) 
nuiic) = H<iUC) 

hlUC) - Hk 

n2UC> = Hi 
faO TO llbO 
1140 CONTINUE 
m3UC)=Hk 
h4 ( IC)=Hi 

1150 CONTINUE 

1040 NC=NC+1 

TO 10^0 
1030 CONTINUE 

6001 rOKMATUhl#4uy»4H — * ibAfcmH »///) . — 

*P1Tl< 10>6U01> (TlILt ( J> • J=lt*S) 

UO 1U60 u=l»Kl 

L=K-J+1 

1060 toPATEUO»6u0<>) Lt AAl i J> *L'bBl ( U) 

6002 FOKPaKIh » lbX*6HALPHA 4 » 12#4H) = *U21.16»20X» 

1 bH PtT A t f I*if 4H) = » 021 • 16 ) 

6u0b h PKP aT 4 /// ) 

rtPITt < 10#6u0b) 

L)0 1U61 u=i»Kl 
L=K~U+1 

luM *P1TL( 10»6 u0a> L» AA2i J) f L*bR2 ( J) 

6004 EOKPATdh » lbX t 6hGaMma ( » 12 > 4H ) = *U21.16r20Xt 
1 oHUELTa < , I2*4H> = t 1)21 * 16) 

600.6 hOKPATUtt. *lbXf5HThETA*17X*5hwbAHl»4bXf ‘SHHbAH?) 

wP lTt.< IOf 6u0b> 

*PiTt< 10»6oQb) 

6U03 FOKMATUm *1 3X t ElO .4 f *sX * 1H { , E 10 . 4 r 3H , , £1 0 . 4 f 1H ) # 1 Ox r 
1 lbX* 1H( »EiO*4 f 3H t »El0.4rlh)) 

UO 1070 1=1 » IC 

wPlTt( 10f 6003) _ 

1070 CONTINUE 

CALL OMPbUP 
UO 11 RU 1 = l»b 
116U bCUXU) = llTLtU) 

CALL OUIinML (—1*— 3«0*b*0»—3*Of3*0»lH t PCDX » BCUV » I C * H 1 # H£ ) 

CALL OUInMH 0»-*3*0»3*P r -3f 0r3*0» 1H » PCDX * bCDY t IC »H3|_H4_L 

UO 1140 l=lrlC 
h?(I)=-H*:(i) 

1140 H41I) = -H4tl) 

CALL oUIkMl ( 0#-3.U»3.0*-3*0,3.0> 1H *PCnx»bCDY* lC'HdH2>. 
CALL QUlKML( 0 r-3* U » 3* 0 * -3 . 0 t 3* 0 » 1H t PCDX * bCDY * IC * H3 9 H4 ) 

hFT UKN _ .2 

ENU 



TABLE BII. - LINE PRINTOUTS OF BASHFORTH -MOULTON 
PREDICTOR-CORRECTOR PAIRS 
(a) k = 1, d = 1 


BASHFORTH-MOULTON K=lr D=1 


ALPHA { 

1) 

= .1000000000000000+01 

BETA ( 

1) □ 

,0000000000000000 

ALPHA ( 

0) 

= -.1000000000000000+01 

BETA ( 

0) = 

,1000000000000000+01 


GAMMA ( 

1) 

= ,1000000000000000+01 

DELTA ( 

1) = 

•1000000000000000+01 

GAMMA ( 

0) 

= -.1000000000000000+01 

DELTA ( 

0) = 

.0000000000000000 


THETA 

HBAR1 

HBAR2 

,0000 

( .0000 t 

-.0000 ) 

( -.1000+01 

-.0000 ) 

,3491-01 

( .6055-03 * 

.3486-01) 

( -.1001+01 

-.3486-01) 

,6981-01 

( ,2378-02 » 

,6943-01) 

( -.1002+01 

-.6943-01) 

.1047+00 

( ,5197-02 , 

,1035+00) 

< -.1005+01 

-.1035+00) 

,1396-00 

( ,8886-02 » 

.1367-00) 

( -.1009+01 

-.1367-00) 

,1745-00 

( .1325-01 » 

,1692-00) 

( -.1013+01 

-.1692-00) 

,2094-00 

( ,1808-01 r 

,2007-00) 

( -.1018+01 

-.2007-00) 

.2443-00 

( ,2321-01 r 

.2312-00) 

< -.1023+01 

-.2312-00) 

.2793-00 

( ,2846-01 r 

.2608-00) 

( -.1028+01 

-.2608-00) 

,3142-00 

( .3372-01 r 

.2895-00) 

( -.1034+01 

-.2895-00) 

.3491-00 

( .3889-01 r 

.3173-00) 

( -.1039+01 

-.3173-00) 

.3840-00 

( .4387-01 f 

.3444-00) 

< -.1044+01 

-.3444-00) 

,4189-00 

< .4860-01 » 

.3707-00) 

( -.1049+01 

-.3707-00) 

,4538-00 

( .5305-01 # 

.3963-00) 

( -.1053+01 

-.3963-00) 

.4887-00 

< .5717-01 r 

.4213-00) 

( -.1057+01 

-.4213-00) 

.5236-00 

( .6094-01 * 

,4457-00) 

( -.1061+01 

-.4457-00) 

.5585-00 

( ,6434-01 , 

,4695-00) 

( -.1064+01 

-.4695-00) 

,5934-00 

( .6736-01 r 

,4928-00) 

< -.1067+01 

-.4928-00) 

.6283-00 

( ,6998-01 f 

.5156-00) 

( -.1070+01 

-.5156-00) 

,6632-00 

( ,7221-01 r 

.5380-00) 

( -.1072+01 

-.5380-00) 

.6981-00 

( .7403-01 * 

.5599-00) 

( -.1074+01 

-.5599-00) 

,7330-00 

( ,7546-01 > 

.5814-00) 

( -.1075+01 

-.5814-00) 

.7679-00 

( .7649-01 r 

,6025-00) 

( -.1076+01 

-.6025-00) 

.8029-00 

( .7712-01 » 

.6232-00) 

( -.1077+01 

-.6232-00) 

,8378-00 

( ,7735-01 * 

,6436-00) 

( -.1077+01 

-.6436-00) 

,8727-00 

( ,7719-01 , 

.6636-00) 

( -.1077+01 

-.6636-00) 

.9076-00 

( .7665-01 f 

.6833-00) 

t -.1077+01 

-.6833-00) 

.9425-00 

( .7572-01 * 

.7026-00) 

( -.1076+01 

-.7026-00) 

.9774-00 

( .7441-01 f 

,7216-00) 

{ -.1074+01 

-.7216-00) 

,1012+01 

< .7274-01 » 

.7403-00) 

< -.1073+01 

-.7403-00) 

,1047+01 

( ,7070-01 » 

.7587-00) 

( -.1071+01 

-.7587-00) 

,1082+01 

( .6829-01 r 

.7768-00) 

( -.1068+01 

-.7768-00) 

,1117+01 

< .6553-01 * 

.7946-00) 

( -.1066+01 

-.7946-00) 

.1152+01 

{ .6243-01 r 

.8121-00) 

( -.1062+01 

-.8121-00) 

,1187+01 

( .5897-01 r 

.8294-00) 

( -.1059+01 

-.8294-00) 

.1222+01 

( .5518-01 f 

,8463-00) 

( -.1055+01 

-.8463-00) 

.1257+01 

< .5106-01 > 

.8629-00) 

( -.1051+01 

-.8629-00) 

.1292+01 

( .4661-01 » 

.8793-00) 

( -.1047+01 

-.8793-00) 

,1326+01 

( .4185-01 * 

,8954-00) 

< -.1042+01 

-.8954-00) 

.1361+01 

( .3676-01 t 

,9112-00) 

( -.1037+01 

-.9112-00) 

.1396+01 

( ,3137-01 , 

,9267-00) 

( -.1031+01 

-.9267-00) 



TABLE BD. - LINE PRINTOUTS OF BASHFORTH -MOULTON 


PREDICTOR -CORRECTOR PAIRS - Continued 
(a) k = 1, d = 1 


1431+01 

( .2568-01 9 

.9419-00) 

( -.1026+01 

9 

-.9419-00) 

1466+01 

( ,1969-01 f 

,9568-00) 

( -.1020+01 

9 

-.9568-00) 

1501+01 

( .1341-01 9 

.9715-00) 

( -.1013+01 

9 

-.9715-00) 

1536+01 

( .6844-02 f 

,9859-00) 

( -.1007+01 

9 

-.9859-00) 

1571+01 

( .0000 9 

.1000+01) 

( -.1000+01 

9 

-.1000+01) 

1606+01 

( -.7117-02 r 

.1014+01) 

( -.9929-00 

9 

-.1014+01) 

1641+01 

< -.1450-01 r 

.1027+01) 

( -.9855-00 

9 

-.1027+01) 

1676+01 

( -.2215-01 f 

.1041+01) 

( -.9779-00 

9 

-.1041+01) 

1710+01 

( -.3005-01 * 

.1054+01) 

( -.9700-00 

9 

-.1054+01) 

1745+01 

( -.3820-01 r 

.1066+01) 

( -.9618-00 

9 

-.1066+01) 

1780+01 

( -.4659-01 » 

,1079+01) 

( -.9534-00 

9 

-.1079+01) 

1815+01 

( -.5522-01 f 

,1091+01) 

( -.9448-00 

9 

-.1091+01) 

1850+01 

( -.6408-01 9 

,1103+01) 

( -.9359-00 

9 

-.1103+01) 

1885+01 

( -.7317-01 9 

.1114+01) 

{ -.9268-00 

9 

-.1114+01) 

1920+01 

( -.8248-01 r 

,1125+01) 

< -,9175-00 

9 

-.1125+01) 

1955+01 

{ -.9200-01 r 

,1136+01) 

( -.9080-00 

9 

-.1136+01) 

1990+01 

( -.1017+00 » 

.1147+01) 

( -.8983-00 

9 

-.1147+01) 

2025+01 

( -.1117+00 * 

,1157+01) 

( -.8883-00 

9 

-.1157+01) 

2059+01 

( -.1218+00 * 

.1167+01) 

( -.8782-00 

9 

-.1167+01) 

2094+01 

( -.1321-00 , 

,1177+01) 

( -.8679-00 

9 

-.1177+01) 

2129+01 

< -.1426-00 » 

.1186+01) 

( -.8574-00 

9 

11B6+01 ) 

2164+01 

( -.1533-00 9 

,1196+01) 

< -,8467-00 

9 

-.1196+01) 

2199+01 

( -.1641-00 » 

.1204+01) 

( -.8359-00 

9 

-.1204+01) 

2234+01 

( -.1752-00 , 

.1213+01) 

( -.8248-00 

9 

-.1213+01) 

2269+01 

( -.1863-00 » 

.1221+01) 

( -.8137-00 

9 

-.1221+01) 

2304+01 

( -.1977-00 * 

,1229+01) 

< -.8023-00 

9 

-.1229+01) 

2339+01 

( -.2092-00 f 

.1237+01) 

( -.7908-00 

9 

-.1237+01) 

2374+01 

( -.2208-00 9 

,1244+01) 

< -.7792-00 

9 

-.1244+01) 

2409+01 

( -.2325-00 r 

,1251+01) 

( -.7675-00 

9 

-.1251+01) 

2443+01 

( -.2444-00 9 

,1258+01) 

< -.7556-00 

9 

-,1258+01) 

2478+01 

( -.2564-00 r 

♦1264+01) 

( -.7436-00 

9 

-.1264+01) 

2513+01 

( -.2686-00 9 

,1270+01) 

< -.7314-00 

9 

-.1270+01) 

2548+01 

( -.2808-00 9 

.1276+01) 

( -.7192-00 

9 

-.1276+01) 

2583+01 

( -.2932-00 9 

,1281+01) 

( -.7068-00 

9 

-.1281+01) 

2618+01 

( -.3056-00 9 

.1286+01) 

( -.6944-00 

9 

-.1286+01) 

2653+01 

( -.3181-00 9 

,1291+01) 

( -.6819-00 

9 

-.1291+01) 

2688+01 

{ -.3308-00 9 

.1295+01) 

( -.6692-00 

9 

-.1295+01) 

2723+01 

( -.3435-00 . 

,1299+01) 

( -.6565-00 

9 

-.1299+01) 

2758+01 

( -.3563-00 9 

,1303+01) 

( -.6437-00 

9 

-.1303+01) 

2793+01 

( -.3691-00 9 

,1306+01) 

( -.6309-00 

9 

-.1306+01) 

2827+01 

( -.3820-00 9 

,1310+01) 

( -.6180-00 

9 

-.1310+01) 

2862+01 

( -.3950-00 9 

.1312+01) 

( -.6050-00 

9 

-.1312+01) 

,2897+01 

( -.4080-00 » 

,1315+01) 

( -.5920-00 

9 

-.1315+01) 

,2932+01 

( -.4211-00 ♦ 

,1317+01) 

< -.5789-00 

9 

-.1317+01) 

,2967+01 

( -.4342-00 r 

.1319+01) 

( -.5658-00 

9 

-.1319+01) 

,3002+01 

( -.4473-00 » 

.1320+01) 

( -.5527-00 

9 

-.1320+01) 

,3037+01 

( -.4604-00 9 

.1321+01) 

( -.5396-00 

9 

-.1321+01) 

3072+01 

( -.4736-00 9 

.1322+01) 

C -.5264-00 

9 

-.1322+01) 

,3107+01 

( -.4868-00 9 

.1323+01) 

( -.5132-00 

9 

-.1323+01) 

,3142+01 

{ -.5000-00 9 

.1323+01) 

( -.5000-00 

9 

-.1323+01) 


TABLE Bn. - LINE PRINTOUTS OF BASHFORTH -MOULTON 


ALPHA ( 2) = 
ALPHA ( 1) = 
ALPHA ( 0) = 


GAMMA ( 2) = 
GAMMA { 1) = 
GAMMA ( 0) = 


THETA 
• 0000 
,3491-01 
,6981-01 
• 1047*00 
•1396-00 
•1745-00 
•2094-00 
•2443-00 
•2793-00 
•3142-00 
•3491-00 
•3840-00 
•4189-00 
•4538-00 
•4887-00 
.5236-00 
•5585-00 
.5934-00 
.6283-00 
.6632-00 
•6981-00 
•7330-00 
•7679-00 
.8029-00 
•8378-00 
.8727-00 
.9076-00 
•9425-00 
.9774-00 
.1012+01 
.1047+01 
•1082+01 
.1117+01 
•1152+01 
,1187+01 
•1222+01 
•1257+01 
.1292+01 
.1326+01 


PREDICTOR -CORRECTOR PAIRS - Continued 
(b) k = 2, d = 2 


BASHFORTH-MOULTON K-2 1 D-2 


• 1000000000000000+01 

-• 1000000000000000+01 

•oooooooooooooooo 


BETA ( 2) = .0000000000000000 

BETA ( 1) = ,1500000000000000+01 

BETA ( 0) = -.5000000000000000+00 


. 1000000000000000+01 

-. 1000000000000000+01 

•oooooooooooooooo 


DELTA ( 2) = .5000000000000000+00 
DELTA ( 1) = .5000000000000000+00 
DELTAt 0) = ,0000000000000000 


HBAR1 HBAR2 



.0000 r 

-.0000 ) 

( -.2000+01 


-.0000 ) 


,3710-06 r 

.3491-01) 

( -.1999+01 


-.4248-04) 


.5926-05 , 

.6984-01) 

t -.1995+01 


-.3386-03) 


•2992-04 » 

•1048+00) 

( -.1989+01 


-.1136-02) 


,9421-04 » 

.1398-00) 

( -.1981+01 


-.2670-02) 


.2288-03 r 

.1749-00) 

< -.1971+01 


-.5157-02) 


.4712-03 , 

•2101-00) 

( -.1958+01 


-•8793-02) 


•8656-03 r 

.2453-00) 

( -.1944+01 


-.1374-01) 


,1461-02 , 

•2806-00) 

( -.1928+01 


-.2015-01) 


•2311-02 r 

•3160-00) 

( -.1911+01 


-.2810-01) 


•3468-02 , 

.3513-00) 

< -.1893+01 


-.3767-01) 


.4984-02 » 

.3866-00) 

( -.1874+01 


-.4889-01) 


.6908-02 » 

,4218-00) 

< -.1854+01 


-.6174-01) 


.9278-02 » 

.4568-00) 

( -.1834+01 


-.7618-01) 


.1212-01 r 

.4915-00) 

( -.1813+01 


-.9215-01) 


.1546-01 » 

•5259-00) 

( -.1792+01 


-.1095+00) 


,1928-01 r 

.5598-00) 

( -.1772+01 


-.1282-00) 


.2357-01 » 

.5931-00) 

( -.1751+01 


-•1481-00) 


•2830-01 r 

.6258-00) 

( -.1731+01 


-.1689-00) 


.3342-01 r 

•6578-00) 

< -.1712+01 


-.1906-00) 


,3886-01 » 

.6889-00) 

< -.1692+01 


-.2131-00) 


.4455-01 » 

,7192-00) 

< -.1674+01 


-•2361-00) 


•5041-01 * 

,7485-00) 

< -.1655+01 


-.2596-00) 


•5636-01 t 

.7768-00) 

( -.1638+01 


-.2834-00) 


•6232-01 » 

•8042-00) 

( -.1620+01 


-.3075-00) 


.6819-01 r 

,8305-00) 

< -.1603+01 


-.3317-00) 


.7391-01 f 

•8558-00) 

( -.1586+01 


-.3560-00) 


.7940-01 • 

•8801-00) 

< -.1570+01 


-•3802-00) 


.8458-01 » 

•9035-00) 

( -.1554+01 


-.4044-00) 


.8942-01 9 

.9258-00) 

{ -.1538+01 


-•4285-00) 


.9384-01 » 

,9473-00) 

( -,1522+01 


-.4524-00) 


,9780-01 # 

,9678-00) 

( -.1507+01 


-.4762-00) 


.1013+00 9 

.9875-00) 

( -.1492+01 


-•4997-00) 


.1042+00 » 

•1006+01) 

( -.1476+01 


-•5230-00) 


.1066+00 9 

•1024+01) 

( -.1461+01 


-.5460-00) 


•1084+00 * 

•1042+01) 

< -.1446+01 


-.5688-00) 


•1097+00 9 

•1058+01) 

( -.1431+01 


-.5913-00) 


•1103+00 * 

.1074+01) 

( -.1416+01 


-•6136-00) 


,1103+00 » 

•1090+01) 

< -.1401+01 


-.6355-00) 



TABLE BE. - LINE PRINTOUTS OF BASH FORTH -MOULTON 
PREDICTOR -CORRECTOR PAIRS - Continued 
(b) k = 2, d - 2 


.1361+01 

( .1098+00 , 

.1104+01) 

{ -.1386+01 

f 

-.6571-00) 

.1396+01 

( .1086+00 r 

,1118+01) 

< -.1371+01 

9 

-.6784-00) 

.1431+01 

[ ,1068+00 # 

.1132+01) 

( -.1355+01 

9 

-.6995-00) 

.1466+01 

( ,1044+00 , 

.1145+01) 

( -.1340+01 

9 

-.7202-00) 

.1501+01 

( .1014+00 r 

.1157+01) 

( -.1325+01 

9 

-.7405-00) 

.1536+01 

( ,9779-01 * 

.1169+01) 

( -.1309+01 

9 

-.7606-00) 

.1571+01 

( ,9362-01 » 

.1180+01) 

( -.1294+01 

9 

-.7804-00) 

.1606+01 

( .8888-01 » 

.1191+01) 

{ -.1278+01 

9 

-.7998-00) 

.1641+01 

( .8358-01 9 

,1202+01) 

{ -.1262+01 

9 

-.8189-00) 

.1676+01 

( .7773-01 9 

.1212+01) 

( -.1246+01 

9 

-.8377-00) 

.1710+01 

< .7135-01 » 

,1222+01) 

( -.1230+01 

9 

-.8561-00) 

.1745+01 

( ,6444-01 * 

.1231+01) 

( -.1214+01 

9 

-.8743-00) 

.1780+01 

< .5701-01 » 

.1240+01) 

{ -.1198+01 

9 

-.8921-00) 

.1815+01 

( ,4909-01 t 

,1248+01) 

( -.1181+01 

9 

-.9095-00) 

.1850+01 

( .4068-01 t 

.1257+01) 

( -.1165+01 

9 

-.9267-00) 

.1885+01 

( ,3180-01 t 

,1264+01) 

< -.1148+01 

9 

-.9435-00) 

,1920+01 

( .2246-01 » 

.1272+01) 

( -.1132+01 

9 

-.9600-00) 

,1955+01 

< .1268-01 * 

.1279+01) 

< -.1115+01 

9 

-.9761-00) 

.1990+01 

( .2462-02 » 

,1286+01) 

( -.1098+01 

9 

-.9919-00) 

.2025+01 

< -.8176-02 t 

.1292+01) 

( -.1081+01 

9 

-.1007+01) 

,2059+01 

( -.1922-01 * 

,1298+01) 

( -.1064+01 

9 

-.1023+01) 

.2094+01 

{ -.3066-01 > 

.1304+01) 

{ -.1046+01 

9 

-.1037+01) 

.2129+01 

< -.4248-01 » 

.1309+01) 

( -.1029+01 

9 

-.1052+01) 

.2164+01 

( -.5466-01 t 

.1314+01) 

( -.1011+01 

9 

-.1066+01) 

,2199+01 

( -.6721-01 * 

.1319+01) 

( -.9937-00 

9 

-.1080+01) 

.2234+01 

( -.8010-01 r 

.1323+01) 

( -.9760-00 

9 

-.1093+01) 

.2269+01 

( -.9332-01 * 

.1327+01) 

( -.9582-00 

9 

-.1106+01) 

.2304+01 

( -.1069+00 r 

.1331+01) 

( -.9404-00 

9 

-.1119+01) 

.2339+01 

( -.1207+00 r 

.1335+01) 

( -.9224-00 

9 

-.1132+01) 

,2374+01 

{ -.1348-00 » 

.1338+01) 

( -.9044-00 

9 

-.1144+01) 

.2409+01 

( -.1493-00 r 

.1341+01) 

( -.8863-00 

9 

-.1155+01) 

.2443+01 

( -.1640-00 t 

.1343+01) 

( -.8681-00 

9 

-.1167+01) 

.2478+01 

( -.1789-00 * 

.1345+01) 

{ -.8498-00 

9 

-.1178+01) 

.2513+01 

( -.1942-00 » 

.1347+01) 

< -.8315-00 

9 

-.1189+01) 

.2548+01 

{ -.2096-00 r 

.1348+01) 

( -.8132-00 

9 

-.1199+01) 

.2583+01 

{ -.2253-00 * 

.1350+01) 

C -.7948-00 

9 

-.1209+01) 

.2618+01 

( -.2412-00 » 

.1350+01) 

( -.7764-00 

9 

-.1219+01) 

.2653+01 

< -.2574-00 , 

.1351+01) 

( -.7579-00 

t 

-.1228+01) 

.2688+01 

t -.2737-00 » 

.1351+01) 

( -.7395-00 

9 

-.1237+01) 

.2723+01 

( -.2902-00 » 

.1351+01) 

( -.7210-00 

9 

-.1246+01 ) 

.2758+01 

( -.3069-00 r 

.1350+01) 

( -.7024-00 

* 

-.1254+01) 

.2793+01 

( -.3238-00 t 

,1350+01) 

( -.6839-00 

9 

-.1262+01) 

,2827+01 

( -.3408-00 » 

.1348+01) 

( -.6654-00 

9 

-.1270+01) 

,2862+01 

i -.3580-00 r 

.1347+01) 

( -.6469-00 

9 

-.1277+01) 

,2897+01 

( -.3754-00 r 

.1345+01) 

( -.6284-00 

9 

-. 1284+01 ) 

.2932+01 

< -.3928-00 f 

.1343+01) 

( -.6099-00 

9 

-.1291+01) 

.2967+01 

( -.4104-00 t 

,1340+01) 

( -.5915-00 

9 

-. 1297+01 ) 

,3002+01 

( -.4281-00 * 

.1338+01) 

( -.5731-00 

9 

-.1303+01) 

.3037+01 

( -.4460-00 t 

.1334+01) 

( -.5547-00 

9 

-.1308+01) 

,3072+01 

( -.4639-00 t 

.1331+01) 

( -.5364-00 

9 

-.1313+01) 

.3107+01 

( -.4819-00 * 

.1327+01) 

{ -.5182-00 

9 

-.1318+01) 

,3142+01 

< -.5000-00 , 

.1323+01) 

{ -.5000-00 

9 

-.1323+01) 


TABLE BII. - LINE PRINTOUTS OF BASH FORTH -MOULTON 


PREDICTOR -CORRECTOR PAIRS - Continued 
(c) k = 3, d = 3 


BASHFORTH-MOULTON K=3 » D=3 


ALPHA ( 

3) 

= 

.1000000000000000+01 

BETA ( 

3) 

- 

.0000000000000000 

ALPHA { 

2) 

= 

-.1000000000000000+01 

BETA ( 

2) 

- 

. 1916666666666666+01 

ALPHA ( 

1) 


.0000000000000000 

BETA ( 

1) 

□ 

-.1333333333333332+01 

ALPHA ( 

0) 

= 

.0000000000000000 

beta< 

0) 

= 

,4166666666666666+00 


GAMMA ( 

3) 

□ 

,1000000000000000+01 

DELTA ( 

3) = 

.4166666666666666+00 

GAMMA { 

2) 


-.1000000000000000+01 

DELTA ( 

2) = 

.6666666666666665+00 

GAMMA ( 

1) 

- 

.0000000000000000 

DELTA ( 

1) = 

-.8333333333333331-01 

GAMMA ( 

0) 

= 

.0000000000000000 

DELTA ( 

0) = 

.0000000000000000 


THETA 

HBAR1 

HBAR2 

.0000 

( .0000 r 

-.0000 ) 

{ -.2400+01 > 

-.0000 ) 

.3491-01 

( -.6141-07 , 

.3491-01) 

( -.2400+01 r 

.4247-04) 

.6981-01 

{ -.9610-06 * 

,6981-01) 

< -.2400+01 r 

,3383-03) 

.1047+00 

( -.4684-05 f 

,1047+00) 

( -.2400+01 » 

.1134-02) 

.1396-00 

( -.1401-04 » 

.1396-00) 

( -.2399+01 » 

.2659-02) 

,1745-00 

( -.3174-04 r 

.1746-00) 

( -.2399+01 , 

.5124-02) 

.2094-00 

( -.5964-04 * 

,2095-00) 

( -.2397+01 f 

.8705-02) 

.2443-00 

( -.9719-04 r 

,2445-00) 

( -.2395+01 * 

.1354-01) 

.2793-00 

( -.1400-03 # 

.2795-00) 

( -.2392+01 > 

.1972-01) 

,3142-00 

( -.1783-03 r 

.3147-00) 

( -.2387+01 * 

.2728-01) 

.3491-00 

( -.1949-03 , 

.3499-00) 

( -.2380+01 » 

.3618-01) 

.3840-00 

( -.1633-03 r 

,3853-00) 

( -.2371+01 r 

.4632-01) 

.4189-00 

{ -.4534-04 r 

.4209-00) 

( -.2359+01 » 

.5750-01) 

.4538-00 

{ .2105-03 r 

.4567-00) 

{ -.2345+01 r 

.6946-01) 

.4887-00 

( ,6720-03 * 

,4928-00) 

( -.2327+01 * 

.8185-01) 

.5236-00 

( ,1425-02 * 

,5291-00) 

( -.2306+01 » 

.9422-01) 

.5585-00 

( ,2574-02 » 

,5658-00) 

( -.2281+01 t 

.1061+00) 

.5934-00 

( ,4247-02 , 

.6027-00) 

( -.2253+01 # 

.1169+00) 

,6283-00 

( ,6591-02 » 

.6400-00) 

( -.2222+01 * 

.1261-00) 

.6632-00 

( .9770-02 » 

.6775-00) 

( -.2187+01 > 

.1332-00) 

.6981-00 

( ,1396-01 » 

,7152-00) 

( -.2149+01 * 

,1375-00) 

,7330-00 

< ,1935-01 f 

.7530-00) 

( -.2109+01 r 

.1388-00) 

,7679-00 

( ,2611-01 f 

,7905-00) 

( -.2067+01 r 

.1366-00) 

,8029-00 

( ,3439-01 # 

,8277-00) 

( -.2023+01 * 

,1308-00) 

.8378-00 

( ,4428-01 * 

.8642-00) 

{ -.1979+01 * 

.1214+00) 

.8727-00 

( ,5581-01 » 

,8997-00) 

( -.1935+01 t 

.1084+00) 

.9076-00 

( .6894-01 > 

,9337-00) 

( -.1892+01 * 

.9211-01) 

.9425-00 

( .8352-01 > 

.9661-00) 

( -.1850+01 r 

.7292-01) 

,9774-00 

( ,9932-01 r 

.9963-00) 

( -.1809+01 # 

.5123-01) 

.1012+01 

( ,1161+00 * 

.1024+01) 

( -.1769+01 » 

,2748-01) 

,1047+01 

( .1334-00 » 

.1050+01) 

( -.1732+01 t 

.2118-02) 

.1082+01 

( ,1509-00 , 

.1073+01) 

{ -.1696+01 f 

-.2445-01) 

.1117+01 

< .1683-00 f 

,1093+01) 

( -.1662+01 r 

-.5186-01) 

.1152+01 

< ,1852-00 f 

.1111+01) 

( -.1629+01 * 

-.7980-01) 

.1187+01 

( ,2014-00 » 

.1126+01) 

( -.1599+01 » 

-.1080+00) 

.1222+01 

( .2165-00 * 

.1139+01) 

( -.1569+01 r 

-.1364-00) 

.1257+01 

( .2304-00 , 

,1150+01) 

( -.1541+01 r 

-.1646-00) 


TABLE BII. - LINE PRINTOUTS OF BASHFORTH -MOULTON 


PREDICTOR-CORRECTOR PAIRS - Continued 
(c) k = 3, d = 3 


1292+01 

( .2429-00 f 

.1159+01) 

( -.1514+01 

9 

-.1927-00) 

1326+01 

( ,2539-00 » 

.1167+01) 

( -.1487+01 

9 

-.2204-00) 

1361+01 

( ,2633-00 t 

.1173+01) 

( -.1462+01 

9 

-.2478-00) 

1396+01 

( ,2712-00 # 

,1178+01) 

( -.1437+01 

9 

-.2749-00) 

1431+01 

( .2775-00 # 

,1102+01) 

< -.1413+01 

9 

-.3015-00) 

1466+01 

< ,2823-00 # 

,1185+01) 

< -.1390+01 

9 

-.3276-00) 

1501+01 

( .2855-00 9 

,1187+01) 

< -.1367+01 

9 

-.3533-00) 

1536+01 

( ,2872-00 r 

.1188+01) 

( -.1344+01 

9 

-.3784-00) 

1571+01 

C .2876-00 t 

.1189+01) 

< -.1322+01 

9 

-.4031-00) 

1606+01 

( ,2865-00 9 

.1190+01) 

( -.1300+01 

9 

-.4273-00) 

1641+01 

< ,2842-00 * 

,1190+01) 

{ -.1278+01 

9 

-.4511-00) 

1676+01 

( .2807-00 r 

.1190+01) 

< -.1257+01 

9 

-.4743-00) 

1710+01 

( ,2760-00 r 

.1190+01) 

( -.1236+01 

9 

-.4971-00) 

1745+01 

( ,2703-00 r 

.1190+01) 

( -.1215+01 

9 

-.5194-00) 

1780+01 

( ,2635-00 9 

.1189+01) 

( -.1194+01 

9 

-.5412-00) 

1815+01 

( ,2558-00 t 

,1189+01) 

( -.1173+01 

9 

-.5626-00) 

1850+01 

( ,2472-00 * 

.1188+01) 

( -.1152+01 

9 

-.5835-00) 

1885+01 

< ,2377-00 » 

,1187+01) 

( -.1131+01 

9 

-.6040-00) 

1920+01 

( .2275-00 * 

.1186+01) 

{ -.1110+01 

9 

-.6240-00) 

1955+01 

( .2165-00 * 

,1185+01) 

( -.1090+01 

9 

-.6436-00) 

1990+01 

( ,2049-00 9 

.1185+01) 

( -.1069+01 

9 

-.6628-00) 

2025+01 

( ,1926-00 9 

.1184+01) 

( -.1049+01 

9 

-.6815-00) 

2059+01 

( .1797-00 , 

,1182+01) 

( -.1028+01 

9 

-.6998-00) 

2094+01 

( ,1663-00 9 

,1181+01) 

( -.1007+01 

9 

-.7177-00) 

2129+01 

( .1523-00 # 

,1180+01) 

( -.9866-00 

9 

-.7351-00) 

2164+01 

< ,1379-00 # 

,1179+01) 

( -.9659-00 

9 

-.7522-00) 

2199+01 

( ,1229+00 r 

.1178+01) 

( -.9453-00 

9 

-.7689-00) 

2234+01 

( .1076+00 » 

.1176+01) 

{ -.9246-00 

9 

-.7851-00) 

2269+01 

( .9183-01 * 

.1175+01) 

( -.9039-00 

9 

-.8010-00) 

2304+01 

( .7568-01 9 

.1173+01) 

{ -.8831-00 

9 

-.8164-00) 

2339+01 

( ,5918-01 r 

,1171+01) 

( -.8623-00 

9 

-.8315-00) 

2374+01 

( .4234-01 r 

.1169+01) 

( -.6415-00 

9 

-.8462-00) 

2409+01 

( .2519-01 » 

.1167+01) 

( -.8207-00 

9 

-.8605-00) 

2443+01 

( .7740-02 * 

.1165+01) 

( -.7999-00 

9 

-.8744-00) 

2478+01 

( -.9988-02 r 

.1163+01) 

( -.7790-00 

9 

-.8880-00) 

2513+01 

( -.2798-01 » 

.1160+01) 

( -.7581-00 

9 

-.9012-00) 

2548+01 

( -.4621-01 » 

.1158+01) 

( -.7372-00 

9 

-.9140-00) 

2583+01 

( -.6468-01 r 

.1155+01) 

( -.7163-00 

9 

-.9264-00) 

2618+01 

( -.8336-01 9 

.1152+01) 

( -.6953-00 

9 

-.9385-00) 

2653+01 

( -.1022+00 r 

.1149+01) 

( -.6744-00 

9 

-.9502-00) 

2688+01 

< -.1213+00 r 

.1145+01) 

( -.6534-00 

9 

-.9616-00) 

2723+01 

i -.1406-00 » 

,1142+01) 

( -.6324-00 

9 

-.9726-00) 

2758+01 

( -.1600-00 9 

.1138+01) 

{ -.6114-00 

9 

-.9833-00) 

2793+01 

( -.1796-00 f 

,1134+01) 

( -.5905-00 

9 

-.9936-00) 

2827+01 

( -.1993-00 9 

.1129+01) 

( -.5695-00 

9 

-.1004+01) 

2862+01 

( -.2192-00 9 

,1125+01) 

( -.5485-00 

9 

-.1013+01) 

2897+01 

( -,2391-00 # 

.1120+01) 

( -.5276-00 

9 

-.1022+01) 

2932+01 

< -.2592-00 > 

,1115+01) 

< -.5066-00 

9 

-.1031+01) 

2967+01 

( -.2794-00 9 

,1110+01) 

< -.4857-00 

9 

-.1040+01) 

3002+01 

( -.2998-00 9 

,1104+01) 

( -.4649-00 

9 

-.1048+01) 

3037+01 

( -,3202-00 9 

,1098+01) 

( -.4440-00 

9 

-.1056+01) 

3072+01 

( -.3406-00 9 

,1092+01) 

( -.4232-00 

9 

-.1064+01) 

31 07,+ 01 

( -.3612-00 9 

,1085+01) 

{ -.4025-00 

9 

-.1071+01) 

3142+01 

( -.3818-00 9 

.1079+01) 

( -.3818-00 

9 

-.1079+01) 


TABLE BB. - LINE PRINTOUTS OF BASHFORTH -MOULTON 
PREDICTOR -CORRECTOR PAIRS - Continued 
(d) k = 4, .d = 4 


BASHFORTH-MOULTON D=4 


ALPHA ( 

4) 

z 

,1000000000000000+01 

BETA( 

4) 

r 

,0000000000000000 

ALPHA ( 

3) 

z 

*.1000000000000000+01 

BETA( 

3) 

z 

.2291666666666666+01 

ALPHA ( 

2) 

- 

.0000000000000000 

BETA< 

2) 

= 

-.2458333333333332+01 

ALPHA < 

1) 

z 

.0000000000000000 

BETAt 

1) 

= 

, 1541666666666666+0 1 

ALPHA ( 

0) 

- 

,0000000000000000 

BETA( 

0) 

= 

-,3750000000000000+00 


GAMMA < 

4) 

z 

.1000000000000000+01 

DELTA ( 

4) 

z 

,3750000000000000+00 

GAMMA ( 

3) 

□ 

*, 1000000000000000 + 01 

DELTA ( 

3) 

= 

. 7916666666666665+00 

GAMMA ( 

2) 

z 

,0000000000000000 

OELTAf 

Z) 

z 

-,2083333333333 332+00 

GAMMA ( 

1) 

z 

.0000000000000000 

DELTA ( 

1) 

z 

.4166666666666666-01 

GAMMA ( 

0) 

= 

,0000000000000000 

DELTA ( 

0) 

= 

,0000000000000000 


THETA 

HBARi 

HBAR2 

0000 

< ,0000 » 

-.0000 ) 

( -.2667+01 f 

-.0000 ) 

3491-01 

( -.2915-09 r 

,3491-01) 

( -.2667+01 t 

.9708-07) 

6981-01 

( -.1854-07 r 

.6981-01) 

( -.2667+01 , 

.3098-05) 

1047+00 

( -,2089-06 r 

.1047+00) 

{ -.2667+01 f 

.2341-04) 

1396-00 

( -.1157-05 r 

.1396-00) 

( -.2667+01 » 

.9800-04) 

1745-00 

( -.4327-05 , 

,1745-00) 

( -.2668+01 , 

.2966-03) 

2094-00 

( -,1262-04 r 

,2094-00) 

( -.2668+01 , 

.7306-03) 

2443-00 

( -.3092-04 , 

,2443-00) 

l -.2670+01 t 

.1561-02) 

2793-00 

( -.6665-04 r 

.2792-00) 

( -.2672+01 » 

.3003-02) 

3142-00 

( -.1300-03 , 

.3142-00) 

( -.2675+01 * 

.5333-02) 

3491-00 

( -.2343-03 , 

.3491-00) 

{ -.2678+01 f 

.8888-02) 

3840-00 

( -.3952-03 r 

.3841-00) 

( -.2683+01 f 

,1407-01) 

4189-00 

( -.6307-03 * 

.4191-00) 

( -.2688+01 * 

.2133-01) 

4538-00 

( -.9594-03 r 

,4542-00) 

( -.2694+01 , 

.3120-01) 

4887-00 

( -.1399-02 r 

.4895-00) 

( -.2700+01 t 

.4422-01) 

5236-00 

( -.1964-02 f 

,5249-00) 

( -.2706+01 » 

.6101-01) 

5585-00 

( -.2663-02 » 

.5607-00) 

( -.2712+01 , 

,8215-01) 

,5934-00 

( -.3493-02 » 

.5969-00) 

( -.2716+01 r 

.1082+00) 

,6283-00 

( -.4437-02 r 

.6336-00) 

( -.2718+01 , 

.1397-00) 

,6632-00 

( -.5454-02 t 

,6709-00) 

( -.2716+01 , 

.1769-00) 

,6981-00 

( -.6474-02 t 

,7091-00) 

( -.2709+01 r 

,2200-00) 

,7330-00 

( -.7382-02 , 

,7484-00) 

( -.2696+01 r 

.2687-00) 

7679-00 

{ -.8000-02 r 

,7890-00) 

( -.2673+01 f 

.3221-00) 

8029-00 

< -.8068-02 , 

.8312-00) 

( -.2640+01 r 

.3790-00) 

8378-00 

{ -.7204-02 , 

.8753-00) 

( -.2595+01 > 

.4371-00) 

8727-00 

( -.4866-02 , 

.9217-00) 

( -.2536+01 , 

.4936-00) 

9076-00 

( -.2893-03 , 

.9704-00) 

< -.2463+01 » 

.5449-00) 

9425-00 

( .7564-02 , 

.1022+01) 

( -.2378+01 , 

.5871-00) 

9774-00 

( ,2004-01 » 

.1075+01) 

( -.2281+01 » 

.6166-00) 

1012+01 

< ,3871-01 # 

.1131+01) 

( -.2178+01 r 

,6303-00) 

1047+01 

( .6515-01 » 

.1186+01) 

( -.2073+01 # 

,6271-00) 

1082+01 

( .1004+00 » 

,1239+01) 

( -.1970+01 , 

.6076-00) 

1117+01 

( .1443-00 t 

.1286+01) 

( -.1876+01 , 

.5744-00) 

1152+01 

( .1954-00 » 

.1325+01) 

( -.1790+01 , 

.5315-00) 

1187+01 

( ,2506-00 r 

.1352+01) 

( -.1716+01 , 

.4828-00) 



TABLE BII. - LINE PRINTOUTS OF BASHFORTH -MOULTON 


, 1222+01 

,1257+01 

,1292+01 

,1326+01 

,1361+01 

,1396+01 

,1431+01 

.1466+01 

.1501+01 

,1536+01 

.1571+01 

,1606+01 

,1641+01 

,1676+01 

,1710+01 

,1745+01 

,1780+01 

.1815+01 

.1850+01 

.1885+01 

.1920+01 

,1955+01 

.1990+01 

,2025+01 

,2059+01 

,2094+01 

,2129+01 

.2164+01 

.2199+01 

,2234+01 

,2269+01 

.2304+01 

,2339+01 

.2374+01 

,2409+01 

,2443+01 

.2478+01 

.2513+01 

,2548+01 

.2583+01 

.2618+01 

.2653+01 

.2688+01 

,2723+01 

.2758+01 

.2793+01 

.2827+01 

.2862+01 

.2897+01 

,2932+01 

,2967+01 

,3002+01 

,3037+01 

.3072+01 

.3107+01 

,3142+01 


PREDICTOR-CORRECTOR PAIRS - Concluded 
(d) k = 4, d = 4 



,3068-00 » 

.1369+01) 

< -.1651+01 

9 

.4313-00) 


,3608-00 f 

,1374+01) 

< -.1594+01 

9 

.3795-00) 
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(a) k = 1, d = 1 

Figure Bl. - FORTRAN IV subroutine STBL2 plots of 
Bashforth- Moulton predictor- corrector pairs. 


50 






-3-2-10123 

(c) k = 3, d = 3 
Figure Bl. - Continued. 
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(d) k = 4, d = 4 
Figure Bl. - Concluded. 
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